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Abstract 

We present detailed discussions on the stochastic Hamihonians for non-critical string field 
theories on the basis of matrix models. Beginning from the simplest c = case, we derive 
the explicit forms of the Hamiltonians for the higher critical case /c = 3 (which corresponds 
to c = —22/5) and for the case c = 1/2, directly from the double-scaled matrix models. In 
particular, for the two-matrix case, we do not put any restrictions on the spin configurations 
of the string fields. The properties of the resulting infinite algebras of Schwinger-Dyson 
operators associated with the Hamiltonians and the derivation of the Virasoro and al- 
gebras therefrom are also investigated. Our results suggest certain universal structure of 
the stochastic Hamiltonians, which might be useful for an attempt towards a background 
independent string field theory. 
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1. Introduction 



A common idea towards a non-perturbative formulation of string theory is to start from the 
concept of string fields. Just as the ordinary local fields describe the motion and interaction 
of particles in terms of creation and annihilation operators, we can construct string field 
theories by appropriately slicing the world-sheets of strings and introducing the field oper- 
ators to create and annihilate the strings. Clearly, there is continuously infinite amount of 
arbitrariness in choosing slicings. For instance, the light-cone string field theory uses the 
light-like plane in the target spacetime to slice the world-sheet, while the covariant string 
field theories [Q], in general, use different methods of slicing, which are based on the geometry 
of the moduli space of Riemann surfaces. The arbitrariness of slicing may be interpreted as a 
sort of gauge freedom of the theory. At present, however, we have no satisfactory framework 
to formulate such a gauge structure in a systematic and general way. 

Recently, an interesting new way of slicing has been proposed Q , and the corresponding 
string field theories @, |^ have been suggested for the case of noncritical strings with c = 
and the case with minimal conformal matter c = 1 — . In this proposal, the world 

sheets are sliced by using a certain time parameter, which is intrinsically defined on the 
world sheet as a measure of the distances from the boundaries of the world sheets. In 
this paper, we will call the string field theories of this type as "proper time" string field 
theories (PSFTs), in analogy with the familiar proper-time representation of propagators 
in ordinary field theories. In PSFTs, it seems less difficult to incorporate the higher-genus 
(and hopefully, non-perturbative) effect than in the moduli space approach as employed 
in the ordinary covariant string field theories. Remarkably enough, there exists a single 
exact Hamilonian operator which directly characterizes all of the correlation functions in the 
system at once. 

One of the many unsolved problems of present PSFTs, however, is that we do not know 
definite symmetry principles, if any, on the basis of which one can more or less uniquely 
characterize the theories. Thus, most of previous attempts had to rely upon certain guess 
works, and one can only justify the theories by checking the agreement of amplitudes with 
known results obtained from other methods, such as the matrix models. In this situation, 
the observation that the PSFT for c = can be interpreted as the collective field theory 
of matrix models formulated in stochastic quantization seems very useful and suggestive. In 
connection with this, we should recall an attractive idea |^ of relating the renormalization 
group formulation of string-field equations to stochastic quantization. 

Another crucial question of PSFTs for further developing the theory is whether or not 
this method of slicing is meaningful for constructing string field theories for c > 1 and critical 
strings. The simplicity of proposed PSFTs for the case c < 1 is of course due to the simplicity 
of the target spaces. For example, in the case where the target space is the Ising model, one 
can deform the slicings such that the spin configuration on each string field is either all spin 
up or all spin down [g]. If one goes to c > 1, the slicing of this type would, however, be too 
singular to be tractable and one would have to introduce string fields without making any 
restrictions on possible matter configurations. 
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From this view point, it seems important to treat even the cases c < 1 without such 
restrictions and to study the structure of resulting PSFT, since we naturally expect that 
such a formulation should exhibit certain universal structure of the general PSFT which is 
common to PSFT for general critical strings. Since there is no known symmetry principle 
on the basis of which we can derive the theories, it is natural to directly derive such a 
formulation starting from the matrix models. That is what we shall present in this paper. 
Our hope is to get some insight on the nature of the PSFTs by deriving the formalism from 
the matrix models as explicitly as possible. We will follow the suggestion of ref. 0, using 
a slightly different approach, and construct the PSFT Hamiltonians directly by taking the 
double scaling limit of the matrix model Hamiltonians.Q 

In the next section, we will first review our method of deriving the stochastic Hamiltonian 
from the matrix models. We illustrate the method by using a simple quantum mechanical 
model with two degrees of freedom and point out some crucial assumptions required for 
proper-time string field theories. In section 3, we treat the case of one-matrix model and 
derive the Hamiltonians for the cases of = 2 (c = 0) critical point and, as a simplest 
example of higher critical models, k = 3 critical point. In section 4, we discuss the Virasoro 
algebra structure associated with the Hamiltonians. Using the example with k = 3, we will 
clarify how the closed Virasoro algebra is obtained for higher critical cases. In section 5, we 
proceed to discuss the two-matrix model. Technically, this case is much more complicated 
than the case of one-matrix model and requires some new ingredients which have not shown 
up in the case of the one-matrix model. We will exhibit some interesting properties on the 
structure of the stochastic Hamiltonians, which may indeed be regarded as an example of the 
universal structure of the general PSFTs. In section 6, we will discuss the closure property of 
the infinite algebras associated with our Hamiltonians. Then, in section 7, the algebra of 
the two matrix model will be derived starting from the infinite algebra. These two sections 
provide consistency checks for the results of section 5, by deriving the expected properties of 
the two-matrix model from the present formalism. In the final section, we will conclude the 
paper by discussing possible implications of our work and remaining problems. Throughout 
this paper, we had to perform a number of tedious computations for which we could not find 
any appropriate references. Most of such details will be described in the Appendix. 



2. The Hamiltonian of Stochastic Quantization 

In this section, we will briefly introduce our method for deriving the Hamiltonian of PSFT. 
For clarity, we take a simple example of zero-dimensional field theory with two degrees of 
freedom x,y with action S{x,y), 



* For previous works which discuss the possibihty of the PSFT with general matter configurations in the 





continuum formulation, see 



i 
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The idea of stochastic quantization [|1^] can be summarized by introducing the following 
Hamiltonian, 

and the Fokker-Planck equation, 

^^{x,y,T) = -H^{x,y,T). (2.3) 

As is well known, the Fokker-Planck equation describes the statistical evolution of the prob- 
ability distribution function \E'(x,?/,r) for the system described by the stochastic equations 
of motion 

dx dS ,^ 

-T- = -TT + ^i' ^2.4) 
ar ox 

dy dS 

where 771, 772 are Gaussian random noises. In the limit of r — > 00, the solution of the Fokker- 
Planck equation reduces to the ground state 

^ ^ e-^, (2.6) 

satisfying H"^ = under the assumption that e^^^"^ rapidly decreases at infinity, corre- 
sponding to the positivity of the hermitian Laplace operator 

eS/'^He-'/^ = DlD^ + DlD2, (2.7) 

with D, = e-^/2^e^/2,D2 = e-^/2^e^/2. 

The Green function of an arbitrary observable O can be expressed as 

(O) = lim f dxdyO{x,y)'${x,y,T). (2.8) 



When ( p.6|) is assumed to be the unique ground state of the Hamiltonian H, we are entitled 
to suppose that the entire Schwinger-Dyson equation is replaced by a single ground-state 
condition given as 

Mm f dxdyO{x,y)H'^{x,y,T) = 0, (2.9) 
which is, after partial integrations, equivalent to 

dxdy[^e-'^ + -^e-'^]O = 0. (2.10) 
ox ox oy oy 

Using the generating functional with O = q'^^^^-'^v ^ this is rewritten as 

H{J,^)Z[J] = Q, (2.11) 



Z[J] = J dxdye-^^'''y^+-^'''+-^^y, (2.12) 

with 



From the assumption of the uniqueness of the solution for the ground state condition (|2.11 
(that we call Hamilton constraint), we can impose 

TiZ[J] = T2Z[J] = 0. (2.16) 

which are nothing but the general form of the Schwinger-Dyson equation for our system, 

TiZ[J] = y rfxc/y — e-^("'J^)+^i"+"^2^ (2.17) 

T^ZU] = / — e-^("'J')+^^"+-^22/_ (2.18) 
J dy 

Note that the integrability condition is automatically satisfied, 

[Ti,T2] = 5,.,(J, ^) - SyM jj) = 0. (2.19) 

In general, by introducing more general source terms J2i Jifi{x,y), the stochastic Hamil- 
tonian takes the following form with a general set of operators Tj 

H = J2 JiTi. (2.20) 

i 

Then, the assumption of uniqueness of the ground state implies that the partition function 
satisfies 

TiZ[J] = (2.21) 

where 

r.ZW^ AM, |.|)e-E.../. p.22) 

which can be, for appropriate choice of the source terms, expressed as functional differen- 
tial operators in terms of Jj's and are equivalent to the Schwinger-Dyson equations of the 
system. This should be regarded as a fundamental assumption of the method of stochastic 
quantization. We note that in general the algebra of the Schwinger-Dyson operators Tj is 
non-Abelian. 
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HZ^ [ dxdy[-^e-'-^ + |-e-^|-]e^+S. J^f' = Q (2.23) 
J ox ox oil on 



Here, we add an important remark which partly underhes our later discussions. Namely, 
by introducing general source terms, together with this assumption, we can make the for- 
malism background independent. After making a shift of the source function Jj — Jj + 6i^s, 
the Hamiltonian equation ( p.20| ) 

d _Q d d _g d 
dy dy' 

is then recast into the following form 

T^Z + H^Z = 0, (2.24) 

where 

H'Z ^ J dxdy{^ + ^)eS. '^^^ = (^ J.T^)Z. (2.25) 

Thus, under the assumption that the Hamilton equation is equivalent to the Schwinger- 
Dyson equations T^Z = 0, the first term of ( |2.24|) vanishes and the Hamilton equation is 
reduced to H^Z = 0, a form which is formally independent of the starting action S. Here, 
7^°'s are the Schwinger- Dyson operators with the shifted source Jj + 6i^s, or in other words, 
with no bare action. 

Now, to gain a more concrete understanding on the above assumption, let us consider a 
simple example with the bare action 

2 3 

Six, y) = V{x) + V{y) + cxy, V{x) = "^ + 3"^. (2.26) 



The operators Ti , T2 are given as 



d d' d 



Ti = -WT - 9^ - c— + Ju (2.27) 



dJi " dJf dJ< 



2 



d d , , 

97n-2-^7rT + J^- (2-28) 



dJ2 dJi 

On the other hand, from the view point of the Schwinger- Dyson equations, it is easy to check 
that the following set of equations gives a closed recursion equation for (x"): 

^ = J dxdy-^x^'e-^, (2.29) 

= / dxdy-^x'^e-^, (2.30) 
J dy 

= / dxdy-^x'^ye-^, (2.31) 
J ox 

which are obtained from the Ti, T2 constraints by making a power series expansion in Ji, J2 
as 



= 0, (2.32) 

j=o 
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Qn 
Qn+l 

-TiZ[J] 



dJfdJ2 



= 0, (2.33) 

j=o 

= 0, (2.34) 

j=o 



respectively. The closed recursion equation for (x") is obtained by expressing the correlators 
of the form {x""y), (x"?/^) using the last two equations in terms of (x") and by substituting 
the results into the first equation.[] 

However, it is not difficult to see that we cannot derive all of these conditions ( |2.29| )^( |2.31|) 
directly by taking finite order derivatives with respect to Ji, J2 from the single Hamiltonian 
constraint, 

(JiTi + J2T2)Z[J] = 0. (2.35) 



For example, by taking a derivative g^Z^gj^ of ( p.35| ), we obtain the sum of ( |2.33| ) and ( |2.34|) 



but can never obtain ( p.33| ) or ( |2.34| ), separately, by taking any derivatives of finite order. 

Thus, the equivalence of the Hamiltonian constraint (|2.35|) with the Schwinger-Dyson 
equations ( ^.l(j| ) is based on the assumption of the uniqueness of the solution for ( |2.35| ) 
which requires that e~'^/^ rapidly decreases at infinity x, ?/ — > ±00. This uniqueness assump- 
tion amounts to setting certain conditions on the partition function Z[J] which cannot be 
expressed in any finite order of the expansion with respect to the source functions Jj. If the 
appropriate global conditions for the uniqueness are not satisfied, the stochastic Hamiltonian 
would fai]0to give a unique ground state in the limit r 00, and the limit would, in general, 
depend on the choice of the initial state. 

In the case of simple quantum mechanical models, it is relatively easy to identify the 
necessary global conditions. However, in more complex systems such as the double-scaling 
limit of matrix models, it is quite nontrivial to state such conditions, and, in fact, there has 
been no known result replying this question. 

It is clear that the PSFT proposed in ref. ^ is based on the tacit assumptions of similar 
nature as above. Unless the Hamiltonian is able to define a more or less unique ground state 
under the same constraint as for the Schwinger-Dyson equations, the concept of the PSFT 
Hamiltonians would become less significant, since in that case we have to recourse to the 
Schwinger-Dyson equations themselves for the definition of the theory. 

In the following sections, we will discuss the Hamiltonians of PSFT for one- and two- 
matrix models using the above methods, keeping those assumptions in mind. We here men- 
tion that our method can be translated into the language of ref. [^] by making a functional 
transformation 

Z[J] ^ {Z\ = (0| exp(5:^, A)Z[J] (2.36) 



j=o 



t This procedure is essentially the same as the one employed in ref. ||2^ to derive a closed subset of the 
Schwinger-Dyson equations for the two-matrix model. 

•t This is obvious for the case of usual laplacian A — didi. There are an infinite number of polynomial 
solutions for A/ = 0. The uniqueness can be guaranteed under the requirement, say, of normalizability 
condition. 
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The source functions and their derivatives are replaced by the annihilation- and creation- 
string fields 

J^ ^ ^ ^ ^1' (2-37) 
respectively. Then, the correlators are given as 

A d 



{Z\44---€\^) = {X{qj)Z\J]\j^^. (2.38) 



Note that under the above transformation the normal ordering J ' ' ' jj ' ' ' for jj is au- 
tomatically transformed into the one ip"^ ■ ■ ■ ip ■ ■ ■ for ipjip"^ . The Hamilton constraint is thus 
{Z\ : H^ipjip"^) := 0, and the state {Z\ is obtained as {Z\ = lim^_,oo(0|e~^''^'. 



3. PSFTs from One-Matrix Model 

In this section, we derive the stochastic Hamiltonians from the one-matrix model at the 
k = 2 and k = 3 critical points which correspond to the matter central charges c = and 
c = —22/5, respectively. 



3. 1 Stochastic Hamiltonian at c = 

We first treat the case of c = 0. Although this case has already been discussed in ref. 
within the framework of the collective field method, we present some details for the purpose 
of explaining our method which is slightly different from ref. 0. 

The generating functional of the c = one-matrix model is defined by 



Z[J] = - I d'^'M e-^try(M)gj.*^ 
Z J 

Z = /rf^'M e-^t^^(^^), (3.1^ 



V{M) = -M^ - -M^ (3.2) 
2 3 

J-'^ = (3.3) 

$(0 = — tr — - — (3.4) 

is a loop operator and the contour of (^-integral L is chosen to be parallel to the imaginary 
axis such that in the region of the right of L there are no poles of $(C)- The source function 



where 
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J{() can take an arbitrary form as a function on L. The variable ( can be regarded as being 
conjugate to the length of the loop in the sense of Laplace transform. 

We start with the ground state condition of the stochastic Hamiltonian 



-NtTV(M)_ 



d 



where M is expanded by the basis oi N x N hermitian matrices {t"}: 



a=l 



Using the identities 



a 



tiAtiB, 
tiAB, 



(3.5) 



we obtain formulas such as 



V -^^{O = -N—^(C? 



— y tr(M"r)tr(t°- 



d ""^ 1 

-^(C"<f(C)-EC']^tr(M"-i-'^)), 



fc=0 



where the symbol is the so-called combinatorial derivative, defined as 

/(C) - /(C) 



(3.6) 

(3.7) 
(3.8) 

(3.9) 



which appears when two loops merge into a new loop. 

We can then reduce eq. ( p.5|) to a functional differential equation with respect to the 
source J(C), 





n 



nz[j] 

dc 



(3.10) 



S 1 , 2\ 

57(0 ^ 



\i(-9ef-9c 



1 



iV2 Jl 2TTi 



■m 



(3.11) 
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The functional derivative is defined for ( on the contour and acts on the source J{C) 
as a delta function 

^ = 27r.5(C-C')- (3.12) 

when both ( and (' resides on the same contour. Ti is the exact stochastic Hamiltonian 
for the c = PSFT before taking the double-scaling limit. The first term represents the 
splitting process of a loop, while the second represents the merging process of two loops. 
Note that the expression (|3.11|) is normal ordered in the sense that the differential operators 
jjl^ always sit right of J(C)- 

Introducing a lattice spacing a, we now take the continuum limit (the double scaling 
limit) a ^ by defining the scaling variables as 

C = C^{1 + ay), g = g^{l-aH), ^ = a^/^g^^ 



where the critical points are 

C = (v^ + 1) ■ 3^/", g. 



3V4 
6 



The meanings of the variables y, t and g^^ are the Laplace conjugate of loop length, the 
cosmo logical constant, and the string coupling constant, respectively. From the result of the 
disk amplitude [TT|], it can be seen that the contribution of the poles of $(C) accumulates 



to a cut of the interval [-(^3 - 1) ■ S^/^ + 0{a'^), C* - 4 ■ + 0{a^)]. The region 

ReC > C* contains no singularities of $(C)- So, in the scaling limit we can choose as the 
contour L the line K^, — ioo, + ioo], which is mapped to the imaginary axis in y-plane. 

In order to obtain the correct continuum limit, we have to subtract a non- universal part 
from the correlation functions. In the present case, it is required only for the one-point disk 
amplitude. Namely, the connected i^-point function 

w^(Ci,---,a) = ($(Ci)---<f(a))c 

is written as 

mCi,---,a) = a'''/Myi,---,yK) + 0{a^'^^'^/') {K>2) 

where |(C~5'C^) is the non-universal part of the disk amplitude, and w is the universal part 
giving the correct continuum limit. 

Thus we redefine the source J{y) and the functional derivative by 

^ liC-gO + a'/'J-, (3.13) 



SJiO ' 6J{y) 

JiO = C'd-'^Jiy). (3.14) 
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The shift ( ^.1^ ) corresponds to the rescahng of the partition function as Z[J] = exp(/^ ^'-^iOiC' 
gC)/2)Z[J{y)]. Then, H becomes 



J-^oo 27li ""[S.TdlV 4 ^ \ 



'ICO 

fioo 



6J{z) 
where 

T = ^-^t, C=^{l + V3f. 

3(1 + 73)2 ' 12 ^ ^ ^ 

Note that in the merging interaction (namely, the term of the form JJjj) the shift of the 
functional derivative does not contribute because 

d^dc'D,iCX')liz-gz') = 0. (3.15) 

After finite rescalings 

6J{y) 5J[y) 

we have the stochastic Hamiltonian in the continuum theory 

J -100 2711 6J(vV J-ioo2m 

rioo 



SJ{y) 



6J{z) 

Piy) = 3y^-^T, (3.17) 

where the overall factor a^^^C^T^C^^^ was absorbed by a redefinition of the fictitious time. 
This result, which has been already known from ref. essentially coincides with the form 
of the c = non-critical string field theory^ proposed by Ishibashi and Kawai |^ , if one uses 
the Laplace-transformed string fields instead of their loop-length representation. 



3.2 PSFT for a Higher Critical One-Matrix Model 

We next treat a case of higher critical point {k = 3(c = — 22/5)).p] This problem is interesting 
since a naive extension of the c = hamiltonian leads to an apparent contradiction as 



discussed in ref. [15 



8 For a derivation of the c = Hamiltonian directly from dynamical triangulation, see ref. |12| 
^ Extension of the formalism of ref. to higher critical cases has been given in . 
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The k = 3 critical theory is reahzed by the potential of the fourth- degree polynomial 

(3.18) 

where is a real solution of the cubic equation 



V(M) = -i^M^ + ^M^ + — M^), 
^ ' N^ 2 3 20 ^' 



25^73^ - 30(73 + 32 = 0, 



or explicitly 



and 



93 = - 



2^1/3 

25/ 



((8 + 3V6)i/3 + (8-3v^)i/3) 
59l - 2 



92 



(3.19) 



(3.20) 



Note that we do not use the well known even critical potential (of sixth order) at A; = 3 
critical point, in order to avoid a complication caused by the Z2 symmetry.[J] 

Now, let us derive the Hamiltonian for the k = 3 critical case. Considering the generating 
functional (|3.1|) with the potential (|3.18| ), we have the Hamiltonian before taking the scaling 
limit : 



6J{C) 2^'^^\ 



5J{z) 



(3.21) 



where / is the integral over the contour encircling the poles of $(C)- Note that / s/{(') 
corresponds to the insertion of a microscopic loop represented by the operator -^rtrM. In 
the c = case, no such term appears if one uses the third order potential of ( p.2|) , because 
of the formula (PT^). 

Next, we need to identify the non- universal parts of loop operators in the scaling limit 
a —>■ defined by 

N 

(3.22) 



^^1 



C = C(l + aZ/), C* 



-5(?3 + 2 



(3.23) 



II We do not know any previous work discussing the fc = 3 critical point using the quartic potential. For 
a brief explanation of the derivation of the quartic critical potential, see the Appendix A. 
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From the results of the Appendix A, we have 



{m)o = lv'{C) + a'/My) + 0{a'/'), 



( —tiM) 

\n /( 



32 + 25^3 
15 



5 



4 



(3.24) 
(3.25) 



where, in the disk amphtude corresponding to a macroscopic loop, the first term ^V'{C,) is 
the non-universal part, while the second term w{y) denotes the universal one: 



w 



(y) = -leJV - \T"'y + ^T^/^)^/^T^ 



(3.26) 



with T = {2(~^)^t. For the microscopic disk amplitude ( p.25| ), the first two terms represents 
the non- universal part, and the third term is universal. 

Using the above results, we see that the source J{y) and the microscopic loop operator 
Co in the continuum theory should be defined by 



JTln = + ^'''t^T' ^^^^ = «"'''CV~(y), 

dJ{Q 2 SJ{y) 



32 + 25(73 



2nt 5J{C') 15 
The string coupling constant is introduced as 



,4 
'-i 

5 



(3.27) 
(3.28) 



7/2 



(3.29) 



Substituting ( |3.22| ), ( pj.23D and these rescaled expressions into the lattice Hamiltonian 
( |3.21|) , we obtain 



J~ioo zm 



-^C!{y' + lTy') + lc.yOo 



-a 



dy r 



6J{yy 25 

dy' J, , 



J —lOO ZTTi J — loo ZTTt 



6J{z) 



+0{a^), 



(3.30) 



where as in the c = case we chose the line [C* — ioo, + ^oo] as the contour L. Note that 
in the merging term the shift of the derivative jj produces a quadratic term with respect to 
J, but it does not contribute to the leading term of 7-^ as a — 0, because 



N Jl 2m JL 2m 2 



2 2 

« 9st 



—Ay) 

too 2m 



ioo 2m 'lO' 



(3.31) 
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After making the rescalings again, 



5'* ^ 6J{y) 5'* 6J{yy 

Oo--^C:a, gst^\C"9si. (3.32) 
we finally obtain the continuum Hamiltonian 

j-loo zvrz oJ{y) J-ioo Ztti 

+4 /^"^ ^J(y) r ^Jiy)dydy'D^{y,y')J- (3.33) 
J-ioo Zm J-ioo Ztii oJyzj 

~p{y) = 5y' + ^Ty + Oo, (3.34) 

where the overall factor a^/^|C*''^ was absorbed into the fictitious time. 

We emphasize that in this Hamiltonian, the tadpole term — ■^J{y)p{y), is not a 
pure c-number, but contains the y-independent operator Oq. This is in contrast to the c = 
case where the tadpole term consists only of the c-number function. Actually, the "operator" 
part of the tadpole is a misnomer. We should rather call it a kinetic term. The Hamiltonian 
description of the higher critical point requires a kinetic term for infinitesimally small loop, 
in addition to the genuine tadpole corresponding to the c-number part of p. 

In the Appendix B, we will determine the operator Oq using the Schwinger- Dyson equa- 
tions. And in the next section, using this result, we confirm that it is just necessary for 
ensuring the closure of the algebra of the Schwinger- Dyson operators appearing in the Hamil- 
tonian. As is discussed in ref. in trying the extension of the c = Hamiltonian to the 



higher critical case, the integrability condition would not be satisfied if one had naively re- 
placed the p{y) of the c = case with c-number polynomials of higher degree. The authors 
in ref. [|I5| proposed a possible way out, which is, however, different from ours. 

Before concluding this section, we derive the disk amplitude from the PSFT Hamiltonian 
( p. 331 ) and compare with the matrix model result as a consistency check of our result. In the 



sphere approximation, the condition 

6 



-nz[j] 



is reduced to the equation for the one-point function w{y), 



= 0, (3.35) 

j=o 



dy[w{yf - {y' + ^ V + y{Oo)o)] = 0. (3.36) 

By demanding that the cut of w{y) resides only on the real axis, as is required from the 
original definition of the loop operator ( p.4|) , both of w{y) and the expectation value of Oq 
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are uniquely determined as 

w{y) 
{Oo)o 



15 
64' 



^4/3 



(3.37) 
(3.38) 



which coincide with the resuhs (|3.26| ) and (|3.25|) obtained directly without using the Hamil- 
tonian, after taking account the rescalings ( |3.32| ). 



4. Derivation of the Virasoro Constraints 

In this section, we examine the integrability condition of the Schwinger-Dyson operators 
associated with the Hamiltonian equation, 

nz[j] = (4.1) 

for k = 2,3 cases. As a warmup exercise, let us begin from the simplest case of c = 0. 
4-1 The Virasoro algebra at k = 2 

Recalling the discussions of section 2, we rewrite the Hamiltonian Ti. in the form 

n = - r ■^J{y)dyT{y), (4.2) 
d,T(v) = d,T„(v) + p(v), (4.3) 

Thus the Schwinger-Dyson equation associated with the c = Hamiltonian is 

dyT{y)Z[J] = 0. (4.5) 
After a straightforward calculation using the functional derivative 

5J{y) 



5J{y') 

we find a closed algebra for To operators 



2m5{y-y'), (4.6) 



[dyMyi).dyMy2)] = ~gl^dyMdy,-dy,)^—{Toiy,)-Toiy2)), (4.7) 

yi - y2 
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The algebra of dyT{y) is obtained by substituting 

/y 
dy'p{y') 

into ( [4.7|) . Then, using the exphcit expression of p (|3.17| ) we find that the effect of p vanishes 

dy^dy^idy, - dy^)D^{yi,y2) J dyp{y) = 0. (4.8) 
Thus, dyT{y) forms the same closed algebra as dyTo{y) 

[^y,Tiy^),^y,Tiy2)] = -gl^dy^dyAdy, - dy,)^—{T{yr) - T{y,)). (4.9) 

yi - y2 

This agrees with the Laplace-transformed version of the result in ref. p. 
4-2 The Schwinger-Dyson operators at k = 3 

The k = 3 case is less trivial. The only difference from the k = 2{c = 0) case lies in p ( |3.34|) . 
Namely, the c-number part of the tadpole term 5?/^ in p gives a non-vanishing effect 

dyA2{9y^ - dy^)D,{yi, 2/2) J dyp{y) = 2{yi - 7/2), (4.10) 

which is the reason why the naive extension violates the integrability condition. However, p 
is not a pure c-number and contains the operator part Oq. Then, on making the substitution 
as in the c = case, the algebra of dyTQ^y) becomes 

[dy,T{y,),dy,T{y,)] = -gl^dy,dy,{dy, - dy,)^—{T{y^) - T{y,)) 

yi - 2/2 

+2^st(^i -^2) 

+dy,[To{y,),Oo]-dy,[To{y2),Oo]. (4.11) 

Since Oq inserts a microscopic loop, it can be expressed by some local operator, obtained as 
some coefficient of large y expansion of the loop operator jj^^- This is done in the Appendix 
B. The final results are 

a = 2 / ^y^/^-i- (4.12) 
Jc 2m SJ(y) 

= 2 hm r ^e'Vyy^^—, (4.13) 

where the contour C surrounds the negative real axis, and the contour [— ioo,ioo] is under- 
stood to avoid the singularity at the origin to the right. Also, the integral is defined by the 
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analytic continuation using the Beta function and the hmit e +0 must be taken after the 
integration. 

Now in calculating the commutator [To(?/),Oo], we can use ( [4.13|) rather than (^4.12|) , 
because the functional derivative 

-j— = 2m5[y-y) 
6J{y') 

is defined for y, y' on the imaginary axis. As a result of the straightforward calculation 
performing partial integration once, the commutator [To(?/), Oq] becomes 

2 [ <¥_ ,_i/2 1 <^ ^ ^ 



[To{y), Co] = -gL / ^y'-^'^- ^ . (4.14) 

Substituting { WM . and using the formulas in the Appendix B, we see that only the term 

Jc 2tii y — y 

survives after some cancellations. Thus, we obtain 

[T,{y),0,] = -gl^y\ (4.15) 
which makes the algebra ( [4.11| ) for T{y) closed: 

[dyj{y{), dy,T{y,)] = -gl^dyMdy, - (T^/i) - T{y,)). (4.16) 

yi - y2 



We note that this algebra is of course identical with the usual Virasoro algebra ( [B.18| ), 
after taking into account the contribution from the transformation ( p.l9| ). All what we 
have done is merely a check of self consistency. It clarifies, however, how the closure of the 
Schwinger-Dyson operators associated with the Hamiltonian is satisfied for higher critical 
points, owing to the presence of the operator part of p{y) . 



5. c = 1/2 PSFT from Two-Matrix Model 

We now apply our method to the two- matrix model and derive the c = 1/2 PSFT without 
making any restrictions on the spin configurations of the string fields. As we emphasized 
in the Introduction, such a treatment will hopefully reveal certain universal properties of 
the PSFT which are basically independent of the structure of target spaces. This is our 
motivation for performing this analysis in spite of its technical difficulties. 
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5. 1 Stochastic Hamiltonian of the Two-Matrix Model 
The generating functional of the two-matrix model is defined by 
Z\J] = l|rf^^Arf^^5e-V-*, 
Z = j d'^'Ad'^'Be-', 

S = NtT{V{A) + V{B)-cAB), V{A) = ^A^ - ^A^^, 



(5.i: 



°° r dC- da- 



(5.2) 



where the components of the string field $ are defined in terms of the matrix variables as 



N C~A 
$b((t) = — tr 



N"a-B' 
1 / 1 



Ci-Aai-B („.-A(T„ 



B 



fn = 1,2,---). 



The component $^ ('^'b) represents a loop on which only a single spin A (i?) is put, and 
represents a loop on which two domains of A- and 5-spins alternatively appear n times. The 
variables Q, ai (z = 1, 2, . . . , ra) can be regarded as the conjugate variables corresponding 
to lengths on the loop for spinstates A, B, respectively. 

As before, the stochastic Hamiltonian for this system is derived from the identity 

Q = - j d'^'Ad'^'BY. 



a=l 



d 



+ 



d 



d 



dA^ dA^ dB^ OB, 



(5.3) 



By extending the formulas for the one-matrix model, we can arrange eq. (|5.3| ) in the following 
form: 





n 



nz[j] 



S 5 



J ■ K— - J ■ — V — - —J • J ■ A — 



5J , 



6J 6J 



SJ. 



J-T. 



(5.4) 



Here, reflecting the cyclic symmetry of pairs (Ci,crj), the derivative -^j- is defined by 
-FTT? 7 T = 0m,n-{2m) 

X E ^(Cl-C(l))^(^l-^c(l))---^(Cn-Cc(n))'^(t^n -<(«))• (5-5) 

cxyclic permutation 
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Let us explain the meaning of each term, 
i) The first term, "kinetic term", J • (^Kjj^, coming from a part of the product of the 

derivatives of the classical action and the source term, symbolizes the contributions which 
preserve the number of string fields. The first few components of the kinetic operator K are 
as follows. 



da 5 f dC >. 3 



+c 



1 da ^ d ^ ci —C 



(5 5 



f da ( () () \ 

+c ?> - — -.a 



2^ \SJ3{Ci,aXucri,C2:Cr2) 5J3(Ci, (7i, C2, (7, C2, (72), 

. dC / S S 
+cf7r^C\ TT77 7 7 ^ + 



27ri \SJ3{Ci,aiX:Cri,C2:Cr2) 5J3(Ci, c^i, C2, (72, C, cr2)/ 

(5.6) 

where the arrow over d indicates that the derivative acts on the whole functions that follow 
it. 

The structure of the higher components can be inferred from these expressions. Basically, 
each component represents one of the following two processes. The first is the propagation 
of string preserving a spin configuration on a loop, with the loop length being either kept 
fixed or decreased by one- lattice unit. The second is a process in which only a single spin is 
fiipped and the loop length is preserved. For example, let us consider (^Kjj^^ {(i,ai). The 
first and last terms express the former process. Note that, as a special case when a domain 
consists of only one spin, the process can annihilate the domain. The last term represents 
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this. On the other hand, the second and third terms express the latter process, with a single 
spin flip preserving the loop length. 

In the case of [Kjj^^, it is noted that -D„{(Ji,G2)jj^^^ in (Kjj^^ (Ci, C2, (^2) rep- 
resents the following loop when it acts on the partition function, 

s 1/1 1 1 ^ 

-D^(ai,a2)y--— - = -tr 



SJi{Ci,a) N \Ci-Aai-Ba2-BJ' 
Namely, the C2-domain has disappeared in $2(Ci) f^i^ C2, cr2)- 

ii) The second term J ■ (^jj V jj^ , coming from the second derivative of the source term, 

represents processes where a string splits into two. The symbol (^jj V jj^ ^ represents the 
result of splitting of a string with the spin configuration /. The first few components are 

V— (0 = -d^ 



5 J 5J)^^ ' ''5Jb{(tY' 

^ yjj] (Ci,ai,C2,a2) = -2j:(^d,^ + jj^d^ ^ 



+2^c(Ci,C2)tt7^^c(Ci,C2)- ^ 



s s 

+2D^{ai, (^2) YT77 ^Da{crl, 0-2^ 



(5.7) 

We hope the structure of the higher terms containing is self-explanatory from these examples, 
iii) The third term -^J ■ (^J ■ (/\^))) coming from the square of the first derivative of the 
source term, symbolizes processes in which two strings merge into a single string. The symbol 
expresses the result of the two strings with the spin configurations /, J merging into 
a single string : 



6Ja{z) 



^JJbb SJb{s) 
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V / n,A V / A,n 

" 5 

= -dcT.dcMQ^C)jjjT- — ^-T, 

j=l f -'nlS,l, <7l, ,z,aj, - ■ ■ , c,„, an) 

(^^JJ^ (Cl, C^l, ■ ■ ■ , Cn, (^n, (X') = (^/\^^ (^'; Cl, C^l, ■ ■ ■ , Cn, CT„) 
n ^ 

= -d^'Y.^'^,^s{(Tj,a')—— -, 

V / n,l V / l,n 

= j2D,{Q,QDUQX[)j^^, 7 ^ 7 r 

+ T.^s{(Tj,(T[)Dt{aj,(T[)— — — , A/ . A 

j=l OJn+l[{,l, (Ti, ■ ■ ■ , S, (,1, r, (Jj+i, "^nj 

(n = l,2,---) (5.8) 

We again hope that the structure of the generic term is self-explanatory from these examples, 
iv) The last term (the tadpole term), which together with the kinetic term is originated 
from the product of the derivatives of the classical action and the source term, shows the 
processes of the annihilation of a string into nothing: 

d( f da 



JL iTXl JL iTXl 



We here would like to emphasize an important property which characterizes all of the 
above formulas and plays an essential role later in studying the scaling limit. Namely, all 
the processes occur locally with respect to the spin domains. Because of locality, more than 
two domains never be created or annihilated at the same time. As a consequence of this 
rule, only the strings consisting of one domain, $a and $b, can be annihilated into nothing. 
Also, only a single pair of domains can participate in the splitting or merging processes, and 
other domains are left intact. 

5.2 Hamiltonian in Continuum Theory 

Let us now consider the continuum limit of the Hamiltonian ( [5.41 ). As in the one-matrix cases, 
the first task in carrying out this is to identify and to subtract the non-universal parts of the 
disk amplitudes. At this point, a new feature arises. Namely, as we discuss in the Appendix 
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C, the non-universal part of a disk amphtude with a given spin configuration contains, 
in general, the universal parts of the disk amplitudes with simpler spin configurations, in 
addition to the non-universal c-number function which has already appeared in the one- 
matrix model. 

By introducing the connected A;-point correlator in the J = background as 



G 



the generating functional is written as 



Ji, ■ ■ ■ , 4 = A, 5, 1, 2, • • • , 



Z[J] = exp 



J-G« + ij.(J.G(^)) + lj.(J-(J-G(^))) + 



(5.9) 



(5.10) 



Now, the investigation of disk amplitudes in the Appendix C shows that the universal part 
$/ of the operator $/ is obtained by a linear transformation of the following form: 



J 



(5.11) 



where A4ij is a mixing matrix of the universal parts, which is the new feature noticed above, 
and 0/ is the non-universal c-number function. The mixing matrix Aiij is upper-triangular, 
i.e. A4ij = for I < J and is invertible. The first few components of (|5.11|) are 



$a(C) 
<I>b(ct) 
*i(C,cr) 

$2(Cl,cri, (2,0-2) 



'i'B(fT) +0b((t), 

VWci^AiC) + ^Bicr)) + <li(C, ^) + 0i(C, ct) 

-10c(Dc(Cl, C2)$a(C) + D^i^l, ^2)$B(ff)) 
-VWc[D^{Cu C2)($l(C, ^1) + <^'l(C, ^2)) 

+ D^iai, a2)($i(Ci, ^) + '^'i(C2, (t))] 
+*2(Ci, 0-1, C2, cr2) + MCu 0-1, C2, 0-2), 



(5.12) 



where 



MO 



02(Ci,cri, (2,0-2) 



c 2 , 

h -(0- 

3g 3^ 

c(l + 2s - \ 

lOcl 



Wc{C + a)), 



(5.13) 
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and c takes its critical value: c = ^"^y^ , and s = 2 + \fl . Then the connected correlators 
are transformed as 



(1) 



where G-f^... j stands for the universal part of G/^'*... / . 
Thus, by introducing the transformed source 

K 



(5.14) 



the generating functional Z[J] in the continuum theory is obtained by the rescaling Z[J] 
e'^''^Z[J], and takes the form 



Z[j] = exp 



The Hamiltonian acting on Z[j] now becomes 
= HZ[J], 



(5.15) 



(5.16) 



n 



iJM-')-(^KiM^ + 



{JM 



T 



iJM-')-[{Mjj + <p)y{M^.+ 



{JM 



A(^A 
V oJ 



(5.17) 



Next, we arrange eq. ( [5.17[ ) to a simpler form in which the mixing matrix disappears. 
We claim the validity of the following equations: 

-{JM-') ■ [k{M^ + </)) j - {JM-') ■ T 
-{JM-')-(^{M^ + <P)V{M^ + <l>)^ 



{JM-' 
{JM-' 



M- 



6 



V M- 



'6J, 



6J, 

{JM-'] 

{JM-') ■ {{JM-') ■ (A0)) 



AM^ 
6J^ 



\SJ 6JJ 

0, 



(5.19) 

(5.20) 
(5.21) 
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where F is a part of the kinetic operator K, representing only the spin flip process. These 
equations make possible to rewrite the Hamiltonian as 



n = -J- F-^\-J- -^y-^\-—J. J. A-^\\. (5.22 





Justification of the equations ( |5.18D ^( |5.21| ) 



We now present the arguments for establishing the above equations. Firstly, we consider 
the spin flip process in the continuum theory. For the universal parts of the disk amplitudes, 
the loop containing a microscopic domain which consists only of a single flipped spin is 
obtained from the loop containing only macroscopic domains by the following rule: 

-^^W^{C) = -s'i^^a9c#(2)(C,cr), (5.23) 

W^i(Ci;C2,ai) = ."i/^aiy(^)(Ci,a,C2,^i), (5.24) 

' ' ' 5 

where the domain a has been shrunk into the microscopic domain by integration. Here 
we borrow the notations of ref. for the disk amplitudes. For their definitions, see the 
Appendix C. The integral symbol /^o" is used in the sense of the integral of the variable x 
in the continuum theory (a = P*(l + ax)) 

where the contour C encircles around the negative real axis and the singularities of the left 
half plane. stands for the critical value = s(10c)~^/^, and a is a lattice spacing. Note 
that there is a sort of finite renormalization represented by the factor s~^. The derivation 
of the above formulas is given in the Appendix D. 

These relations reflect the property that the spin flip process occurs locally with respect 
to domains; in ( |5.23| ), ( [5. 241 ) only the a domain is concerned and the other domains do not 



change at all. This implies that the relation such as (|5.23| ) and ( |5.24| ) should hold for any 



amplitudes with arbitrary number of handles with generic spin configurations. The complete 
general proof of this important property would be, however, technically formidable, since 
it would require explicit identification of the universal part for general loops with arbitrary 
spin configurations. We can now rewrite the spin flip process in (^K (j^jj + 'i^))^ (C) cis. 




--ddaC + f (C - gC) + J(C - gCYl (5-26) 

c 9g 9 
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where the first term is the universal part, and the others are the non-universal parts which 
are obtained by replacing W{C) with jj^^^ in 
Similarly, we can derive the expressions for 




and 



I^a[(M-L] (Ci,^,Ci,^i) + 02], 
J 2m \ 6J J 2 

f^<l[^Jj) (Ci,ai,C,ai)+02] in (^k(^M^ + <P]] (Ci,o-i). 
Using these results, we arrive at eq. ( |5.18| ) with 




\ 6j) h2m ^^^^ ^ ^'J2m 6 J, 



da i , , _i , - N / '^C ^ ^ 



+ / JB(cT)cs-\-d^)i—C^. 



+ 



r dCi doi c- 1 
/ — — Ji(Ci,cri)cs 
Jl 2m 2m 



7 da 6 
(t> a- 



J 2m 5J2(Ci,^,Ci,^i) 



J 2m 6J2{Ci,auC,(Ti] 
+ •••, (5.27) 

where the ellipsis stands for the terms containing the higher components J„ {n >2). 

We here note that the tadpole term is cancelled with a contribution of the same form 
from the kinetic term. This can be regarded as a consequence of our definition of the disk 
universal part of the single-spin flip amplitude Wi{(): Roughly, the kinetic term contains 
the product of the form J a x one-tpintiip + ^ B), and hence there is always a freedom 

of absorbing the tadpole terms, which are polynomials times the sources Ja (Jb), by appro- 
priately defining the universal parts of the spin-flip disk amplitudes and making the shifts 
for one-spin flip ( c .one-spin flip ) • In thc casc of ouc-matrlx model, the spin-flip process is absent 

so that the tadpole term is directly responsible for determining the disk amplitude. 

Although we do not elaborate further on determining the explicit continuum limits for 
higher components, it is natural, because of the local nature of the spin-flipping process, 
to suppose that the above expression already indicates the generic structure of the kinetic 
term, namely the flipping of a single spin with general spin configurations, the absence of 
the tadpole term and the kinetic term without spin flipping. 
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Secondly, we consider the splitting and merging processes, eqs. (|5.19|) and (|5.2(]| ). In the 
Appendix E, it is explicitly checked that 

for I = A,B,1,2, and that 

for (J, J) = (A A), {B, B), (A, 1), {B, 1), {A, 2), (B, 2), (1, 1). 

As seen from eqs. ( ^.11|) and (|5.12|) , Ai mixes the operators with the configurations which 
are smaller parts of the configuration /, when it acts on the loop operator Eq. ( |5.28| ) 



(eq. (|5.29|) ) means that the mixing is commutative with the splitting (merging) processes 
of loops. This result can again be regarded as a consequence of locality of the splitting and 
merging processes. It is then reasonable to assume that the splitting (merging ) commutes 
with mixing matrix for arbitrary spin configurations. This assertion is nothing but the claims 
( |5.19| ), (|5.20| ). Proving this for completely general case is prohibitively difficult in our present 



technical machinery for treating the double scaling limit. We have to be satisfied with the 
explicit confirmation of this property for several nontrivial cases as given in the Appendix 
E. 

Thirdly, in regard of eq. ( |5.21| ), by direct calculation we can check that 



(A0)/,j = (5.30) 

for (/, J) = {A, A), {B, B), {A, 1), (S, 1), {A, 2), (5, 2), (1, 1). In this case, we can give a gen- 
eral proof of this equation as follows. First we show that (pk niust be polynomial for general 
k. Suppose that this is valid up to some k — 1. Then, from Staudacher's recursion equations, 
as described by ( p.21| ), which relates W^^''^ with amplitudes with lower values of k, we can 
see that the part of the numerator for W^'^''^ consisting only of is a polynomial, because in 
general the combinatorial derivative of a polynomial is also a polynomial. The denominator, 
on the other hand, behaves in the scaling limit as 

Pi - gP^ - cQi - W{Pi) = ~a^{y, + x,) + 0{a'/^). 



Thus, by using the scaled variables Pj = C*(l + aVi), = C*(l + ci^i), the form of (pk can be 
written as 

Polynomial of (?/i , Xi , ■ ■ ■ , , ) 
(Pk = ■ • 

yi + xi 

However, from the cyclic symmetry of (Pj, Qi) in W^'^'^^Pi, Qi, ■ ■ ■ , Pk, Qk), the denominator 
yi + xi must be cancelled with the numerator, and thus 0^ should have the form: 

(pk = Polynomial of {yi,xi, ■ ■ ■ ,yk,Xk), 
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where the polynomial has the same symmetry as W^'^''\ Thus, by induction, the 0^ is a 
polynomial for general k. 

Now from the scaling of the universal part ly^^'^) 

W^^'^Pi, Qi, ■ ■ ■ , Pfc, Qk) = ai-i^w^^^\y^, xi, ■ ■ ■ , y^, x^), 

which is presented in the Appendix C, we expect that the relevant part of 0^ takes the form 

^-{r- III (^-^1) 

Since every component of A0 contains the derivative or the combinatorial derivative, eq. 
( CT) leads to 

{A(j))i^j = for the general components, (5.32) 

as is claimed. 
Continuum Limit 

After these preparations, we are now ready to take the continuum limit of the stochastic 
Hamiltonian ( ^.22|) . From the scaling behaviors of disk amplitudes presented in the Appendix 



C, the scaling of various variables are fixed as follows: 



,2 * 



g = g,{l-a'—T), C = P,{l + ay), fT = P,(l + ax), 



1 

N 



(5Ji(C,a) * SJi{y,x) 

6 



MCu 0^1, C2, cr2) = a~'\J2{yu I/2, X2), 



where the critical values are |]l8l |T^ 



(5.33) 



= VUk^, P, = s(10c)~^/l 

Indeed in the limit a all the universal contributions in the Hamiltonian start with 
0(a^/^), as it should be for the correct continuum limit. After finite rescaling as 

~ 2"" ^ 2^ 2 42 

Ji K Ji, jj^^ P* ^st P* ^St' 
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and absorbing the overall factor a^/^ into the fictitious time, we obtain the continuum Hamil- 
tonian in the form 



6 



6 



6 



""-'■['sj -'■[sj'U -^^^'■['■['Jj^ 



(5.34) 



where the inner product is defined by 

dy 



f-9=[ ^■fA{y)gA{y)+ I ^fB{x)gB{x) 



+ Z] / . n o J o J Xi,---, Vn, Xn)gniyi, Xl, ■ ' ' , Vn, Xn)- 



i J-ioo 1 27ri 27ri 
-i 1=1 

As already emphasized above, only the spin fiip process survives in the kinetic term: 

dx 6 



(5.35) 



;4) - 



c 27X1 6Ji{y,x)' 

dy S 
c 2771 SJi {y, x) ' 
dx 5 

c27ri5J2{yi,x,yi,xi) ' Jc 27ri 5J2{yi,xi,y,xi)\ ' 



+ 



X 



dy 



F-^j {yi,Xi,y2,X2) 



— cs 



+ 



I f dx / 5 ^ 

Jc 27ri \SJ3{yi, x, yi,xi, y2, X2) 
r dy_ ' 
Jc 27ri 



5J3{yi,xi,y2,x,y2,X2) 
5 



SJ3{yi,xi,y,xi,y2,X2) SJ3{yi,xi,y2,X2,y,X2) J \ ' 



(5.36) 



The forms of the splitting ^ V ^ and the merging are the same as in the lattice 
theory, because of the commutativity of the mixing matrix with the processes: The first few 
components of the splitting term are 



-9, 



52 



'SJA{yr' 

5Jb{xY' 
( S 



— V — 
5J 5J , 



{x) = -a^ 



-2 



5JA{y) ^ 5Jb{x) / 5Ji{y,xy 
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\SJ SJJ2 j=i\dJAiyj) SJB{xj) 



+2 (Dy{yi,y2)-~-^ -) I Dy{yi,y2)-~-^ 

V dJi{y,Xi)J V dJi{y,X2), 

+2 \Dx{xi,X2)-~- {Da:{xi,X2) 



5Ji{yi,x)J \ ' 5Ji{y2,x)j' 

(5.37) 



Similarly, the merging processes are given as 



A^) iy;y') = -dydy>D,{y,y')- ^ 



) A,A ^Ja{z) 

A-i) {x]x') = -d^d^>Ds{x,x')- ^ 



j {yi, xi,---,yn, Xn, y') = (^A^ j {y'\ Z/i, • • • , Xn) 



n ^ 

-dy>Y,dy^D,{yj,y') 



j=l ^"^nvZ/l) -^1) ' ' ' J Zj Xj, ■ ■ ■ , ym Xji) 



A^ ) {yi, xi, • • • , y„, Xn, x') = ( A^ ) (x'; yi, xi, • • • , y^, x„) 

n ^ 

= -9,./ d^. D, {Xj , x') -y- , 

j=l ^'Jn\yii Xii' • • 1 yji S, • • • , ym Xn) 



A^ j {yi,Xl, ■■■,yn, Xn] y'l, x[) = (y^JJ^ (^1' ^1' ^1' • • • ' ^n) 

= Y.D^{yj,y[)D^{yj,y[)-j 

j=l OJu+iyyi, Xi, ■ ■ ■ , yj-l, Xj-i, Z, X^, W, Xj, • • ■ , yn, Xn) 

n ^ 

+ Y.^s{xj,x'^)Dt{xj,x\)-j — -, 

j=i o^n+iyyit xi, ■ ■ ■ ,yj, s,yi,t, yj+i, xj^i, • • • , yn, Xn) 

(n = l,2,---) (5.38) 
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We remark that the final Hamiltonian has no tadpole term and no dependence on the 
cosmological constant T. Thus in the two-matrix model case, the cosmological constant 
should be regarded as an integration constant. This contrasts with what one would naively 
expect from the result for the one-matrix cases. In section 7, to check consistency of the 
above results, we will discuss how to obtain a closed set of Schwinger-Dyson equations, 
leading to the W3 constraints, from this Hamiltonian. 



6. Closure of the Schwinger-Dyson Algebra 

Now we proceed to discuss the algebra of the Schwinger-Dyson operators associated with the 
Hamiltonian ( ^.34 ). Comparing with the case of one-matrix models, this requires a much 



more intricate analysis, since there are an infinite number of components for the T operators: 



/too fjni _ rlOO fjq' _ 

^.JAm-dyTi{y)) + / —JB{x){-d.Ti{x)) 
-ioo ZTTZ J—ioo ZTTZ 

+ £/ n ^^^n(2/ba;i, ■ ■ ■ ,2/n,a;„)T„(?/i,a;i, ■ ■ ■ ,2/n,a;„), (6.1) 



(4)), <«'^> 
(4)) j^'- 



- [jj^jj^ {yi,xi,---,yn,Xn) 

~ {^'^ ' {^Jjj) iVi^^ir- ■ ,yn,Xn). (ra = l,2, ■■■) (6.4) 

The Tn operators appearing in the Hamiltonian can be regarded as the symmetrized versions 
of the general Schwinger-Dyson operators T^,T^: 

Ti{yi,xi) = T^{yi;yi,xi) + Tf{yi,xi;xi), 
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T2{yi,Xi,y2,X2) = T^{yi;yi,Xi,y2,X2) +T^iyi,xi,y2]y2,X2) 

+Ti{yi: xi;xi, 1/2, X2) + T^{yi, xi, 1/2, X2; X2), 

(6.5) 

where the semicolon in the argument in the right hand side denotes the point where the 
deformation of a loop occurs in constructing the Schwinger-Dyson equation. For example, 
T^{yi; yi,Xi) {Tf{yi, xi;xi)) represents a deformation of the loop with one pair of domains 
by attaching on it any loops which have at least one A{B)-domedn. The explicit forms of 
the general Schwinger-Dyson operators T^, are 

yij Xi, ■ ■ ■ , yji, Xji) 

+ TTl 7 D^{yi, yn+i)- 



\SJA{yi) SJA{yn+l)J SJn{z,Xi, - ■ ■ ,yn,Xn) 

D^{yi,yj+i)-y- 



— cs 



j—i SJj(^z, Xi, ■ ■ ■ , yj, Xj^ SJji^ji^w, Xj^i, • • • , ym x^) 

I f dx' 5 



c 27ri 6Jn+i{yn+i, x', yi, xi,---, y„, x„) 



/'too (J^y' ^ ^ 
— JA{y')dy>D,{y„+i, y')D^{z, yi)j~- r 
-too Zm dJn{w,Xi,--- ,yn,Xn) 

00 MOO ™ fly', dx'- ~ 



X — 



^Jn+miyii Xii ■ ■ ■ , yj—li Xj—ii Z, Xi, ■ ■ ■ , yn, Xn, W , Xj, ' ' ' ) "^m) 

(6.6) 



-^n {y^t Xli ' ' ' 1 ym Xni 3^n+l) 



+ r D^[xn,yn+i)- 



\5JB(Xn) SJB(Xn+l)J SJn(yi,Xi,---,yn,z) 

""^ 5 S 

^ ] -Pz(-^n+l) Xj) r T / r-Dw{Xj,Xn) 



— CS 



j=i SJj{yi,xi,---,yj,z) 5Jn-j(yj+i,Xj+i,---,yn,w) 

1 / dy' 6 



Ic 27ri SJn+iiyi, xi,---, y„, x„, y', Xn+l) 

/ioo (^2;' ~ 5 
— JB{x')d^'D;,{Xn, x')Dy,{z, Xn+l)-^- 
-ioo 2m dJn{yuXi,---,yn,w) 



°° rioo ™ dy'. dx'. ~ 
5'st / 11 7^'^'^rniyi, Xi,---, y^, X ) 2^ Dz{Xj, Xn+l)Du,{Xj, Xn) 

S 

X— 

SJn+miy'i, x'l, • • • , y'j-i, Xj_i, y'j, z, yi, xi, • • • , yn, w, y'j-^-i, x'j^i, • • • , y'^, x'^, 
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(6.7) 



It should be noted that the Hamiltonian contains the particular set of the above operators 
T^, T^, with partial identification of the variables as yn+i = Vi, Xn = a;„+i, respectively. 

We also remark that, as is expected from our construction of the Hamiltonian, the general 
Schwinger-Dyson operators introduced above are the continuum versions of the Schwinger- 
Dyson operators of the original two-matrix model, appearing in 

= 1 d'^'Ad'^'B T^^^'-^ (C, • • ■)e-^+-'-*. (6.8) 

The explicit forms of the matrix model operators are 



1 ^ d f 1 . 1 



N^^.dAa \Cn+i-A (^-Aai-B Q-Aan-Bj' 

rp{ma.t.)B^^^^ • • • , Cn, CTn] (Tn+l) 

We now compute the algebra of the Schwinger-Dyson operators, as in section 4, by 
changing the contour C to the imaginary axis. For the commutators of the operators with 
the same superscript, the result is found to coincide with that of the corresponding matrix 
model operators, after making the identification: 



They are 



1 . 



{dyn{y),dy,Ti{y')\ = -gl^d,dy,{d, - 5.')-^(ri(z/) - Ti(yO), (6.10) 
[dyT^{y\T,\y'-y\x')] = -^9, |^-^(T,^(y; y', x') + T/^(y'; y, x')) 



i^(T,^(y; y, x') + ?,T,\y'- y', x')) + ^,dy.T^{y'- y' , x') | , (6.11) 
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Similar result is obtained for the operators with the superscript B. 

For commutators between the operators with different superscripts, the situation is not 
so simple, except for the trivial case 

[dyT^iy),d.TE{x)] = 0. (6.12) 

For example, for the commutator [dyT^{y),Tf{y',x';x')], the closed algebra 

[dyT^iy), Tf{y', x'; x')] = gl^dydy> {Tf {y , x'; x') - T^{y\ x'- x')), (6.13) 

which coincides with the result of the matrix model operators, can be obtained only when 
the formula 

f dx 6 

dy DJx,x')DJx,x')^ = (6.14) 

'Jc2m ' ^ "'^ ' '6J2{y,z,y\w) ^ ' 

holds. To check this formula, we need to know a more detailed property of the functional 
derivatives ti^, (n > 2). 

Here we try to justify this formula by using the short-domain expansion of the functional 
derivative operators. It will be useful also for the derivation of the constraints discussed 
in the next section. The short-domain expansions, in general, depend on the choice of a 
background. Since a specific background was already assumed in taking the continuum 
limit of the two-matrix model, use of the short-domain expansion should be allowed here, 
as we have already done in the one-matrix cases in studying the closure of the Schwinger- 
Dyson algebra as a consistency check of the continuum results. The difference between 
the present case and the one-matrix cases is that the tadpole term is absorbed into the 
definition of the spin-flip amplitudes. Because of this peculiarity, the Hamiltonian (|5.34| ) 



itself apparently contains no terms which fix the dimensions of the loop operators. In order 
to fix the dimensions, however, it is sufficient to give the dimension of the Hamiltonian, 
since it then determines the dimensions of the loop operators and the string coupling g^^ 
uniquely. Following our result of a direct derivation of the Hamiltonian from the two-matrix 
model, we assume the dimension of Ti to be \H] = [yY^^- Also, it is noted that a new 
dimensional parameter can enter as the integration constant of the Schwinger-Dyson equation 
dyT^{y)Z[J] = or dxT^{x)Z[J] = 0. We can relate this parameter to the cosmological 
constant with dimension two ([T] = [y]"^), characterizing the c = 1/2 background. See the 
next section, in particular, eq. (|7.20|) , for more details. 

Let us now write down the short-domain expansions. Firstly, for the operators with 



a single domain, jj-, jj-, by considering the dimensions of the loop operators and the 



_5_ 

)Je rr—TT^ 

cosmological constant, we can assume the expansions analogous to the one- matrix case ( |B.2U| ) 
6 



SJAiy) 
5 

6Jb{x) 



W + <^-2/,Ty-"' + E ^nZ/— ^ {y : large), (6.15) 

u 

0^11,^^'^ + a-2/3^a;~2/^ + BuX-""-^ {x : large), (6.16) 
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where a's stand for dimensionless constants, and u runs over positive integers divided by 3. 
The first two terms come from the disk singular terms (thus, af^^ = af^g, 0^2/3 = '^^2/3)' 
and the remaining terms represent both the contributions of the large y or x expansion of 
the cylinder singular parts and those of the local operators. 

Secondly, referring to the explicit form of w^'^\y,x) in ( |C.27|) , we can assume that the 
operator has the following large x expansion: 



S 



SJi{y,x) 



al/^x^/^ + alij^x^l^'y + x-^l\a\i^y' + a'ii/gT) 



(0^7/3?/^ + a%i^Ty'^ + a"}y^B2/^y + a "7/3^5/3) 

+x-'/'{{ai,/,y' + a%^,T)-J— + a%,B,) 
+x-'/'a^%B, ^ 



+x-\a\Or{^){y)+a\B,„) 
oJ 

+a;-2(ai22/Oi(-^)(y) + a%B^/;y + a%B,/s) 
+x-'l'a%0,{jj){y) + 0{x-^/'^), (6.17) 

where Om{jj){y) represents the loop operator with a microscopic m-spin fiipped domain 
being added to the domain y, and the remaining notations are the same as before. Here, 
we assume that the nonpolynomial powers of y all come from the loop operators with the 
macroscopic y domain, and that the fractional powers of T appear only through the S's. 
For every disk amplitude whose explicit form is derived in the Appendix C, we can confirm 
these properties. The above form (|6.17|) is then a consequence of the dimensional analysis. 
Note that we do not have to worry about the ordering of -B's and loop operators, because 
[5.,X]=0due to [^,^]=0. 

Similarly, the large x expansions of I yi x') ' 5j2{y^x y' x) assumed to be, respectively, 

^ ajx + al/^x^^^D^{y,y') ^ 



SJ2{y,x,y',x') SJa{z 
+al/sX^/^D,{y,y" ^ 



6Ji{z, x') 

+aly + a'^y' + a'o^x' + 0{x'^^^), (6.18) 
alx + al/2x'^/^D^{y,y') 



6J2{y,x,y',x) 6Ja{z 

34 



+aly + a'^y' + 0{x-^/'^). (6.19) 

We can fix the normalizations of tlie spin flipped loop operators, referring to the disk 
results in the Appendix D. For example, a5_^ = s. This implies the following relations 

5 _i f dx 6 -1/1 



6 J Jc" 2m Jc 27rz Jc 2Tii 6J2{y, x, y', x') 5JA{y) 

(6.20) 

with d being some constant, which are a consequence of the above expansions. The first 
equation of ( |6.20| ) is the integrated version of eq. ( ^.23] ). Roughly speaking, Om{jj){y) is 
obtained from -j-f- by integrating the domains except y, and the factor is regarded as a 
sort of finite renormalization accompanied with flipping of a single spin. 

Now substituting the expansions ( |6.18| ) and ( |6.19| ), and using the formula such as ( [B.2 



we can directly check the validity of the formula ( |6.14]) , and thus justify the algebra ( |6.13|) . 
From these results, it is expected that the algebra of higher operators also coincides with 
those of the matrix model operators before taking the double-scaling limit. In order to prove 
this for the general case, we need to obtain more precise information with respect to the 
coefficients a's in the expansions of generic loop operators. We would like to emphasize again 
that, since we have started from the matrix model which is already an integral solution of 
the constraints, the closure of the Schwinger-Dyson algebra is guaranteed in our approach. 
The confirmation of the closure of the algebra is useful, however, as a consistency check of 
the continuum limit. Since the structure of the algebra of the Schwinger-Dyson operators is 
essentially determined by the splitting and merging processes of the loops which are closely 
related with those in the Hamiltonian, it is natural to suppose that the algebra is not affected 
by the scaling limit, in view of our discussions in section 5. 

It is straightforward to see that the matrix model operators ( |6.9D form a closed algebra 
such as 



rrp(ma.t.) A rj-i(mat.) A-1 
i-^A ^^A \ 




rT-i(mat.) A 
^A 1 


rrT-,(mat.) B 
[-'■B 5 


r7-i{mat.) B-i 
J-B \ ~ 


rT-i(mat.) B 
^B ) 


rrp{ma.t.)A rp(raa.t.) A-\ 
[^A J-^n J 


rsj 


rpivClsA,.) A 


rrT-,(mat.) B 


7^{mat.)-B] ^ 
n J 


T-i(mat.) B 
^ n 1 


rT-i(mat.) A rp(raa.t.) Ai 
i n ' m J 


rsj 


rT-i(mat.) A 
^ n+m J 


(mat.) B 


7^{mat.)-Bl ^ 
m J 


rT-i(mat.) B 
^ n+m 1 


rrp(vaa.\..) A rp(in-At.) B^ 
[-'-A '-^n J 




T-i{mat.) B 
n 


rry-i(mat.) B 
[ B 


7^ (mat.) A] ^ 
n J 


T-i(mat.) A 
n 1 


rrj-i{mat.) A rj-i{mat.) 

[^A iJ-B \ - 





rm{mat.)j4 T-i{mat.)_B 


rT-i(mat.) A 


. ry-i(mat.) B 

~r n+m 






(n. 


m > 1). 







(6.21) 

The corresponding algebra of the continuum generators 
dyTiiy), T^{y; yi, xi, ■ ■ ■ , x„), 9^T|(x), T^{yi, xi, ■ ■ ■ , ?/„, x„; x) (n = 1, 2, ■ • ■) (6.22) 
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must have the same structure. Half of the generators with the superscript A{B) form a 
subalgebra of the full algebra .0 We note that, as already emphasized in the beginning of this 
section, the generators contained in the Hamiltonian Ti are only the particular symmetrized 
combinations of the general Schwinger-Dyson operators. To ensure the closure of the full 
algebra, we are lead to introduce all of the above general Schwinger-Dyson operators. 



7. Reduction to the W3 Constraints 

In the previous section, we obtained the huge algebra of the Schwinger-Dyson operators. 
Although this algebra has a very complicated structure, we will next demonstrate how the 
W3 constraints, characterizing the c = 1/2 noncritical string, is naturally derived from the 
integrability condition of the first few constraint operators appearing in our Hamiltonian. 
This will provide us yet another consistency check of the preceding results. 
The constraints associated with the Hamiltonian ( |6.1D are 

-dyT^{y)Z[J] = 0, (7.1) 

-d^Ti{x)Z[J] = 0, (7.2) 

T,iy,x)Z[J] = 0, (7.3) 

T2{yi,xuy2,X2)Z[J] = 0, (7.4) 



Considering the closure of this constraint algebra 
[dyTiiy),T,iy',x)] 

[ [y - y ) [y - y r 



+ ^y^y,y (T^(y; y\ x) + Ti^(y'; x) + Tf (y, x; x) + T^{y\ x; x)) | , (7.5) 



which follows from (|6.11|) and ( |6.13|) , we obtain a new combination of the Schwinger-Dyson 
operators, 

(Tfiy; y', x) + T^{y'- y, x) + Tf (y, x; x) + (y', x; x))Z\J\ = 0. (7.6) 
Similarly, from the algebra [dxT§{x),Ti{y,x')], we have 

iT,^iy,x;x')+Tfiy,x' ; x) + T^{y- y, x) + T^{y- y, x'))Z[J] = 0. (7.7) 

** The algebra similar to this subalgebra is presented in ref. M. 
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Let us consider the integrated versions of these: 

dx f dy' 



Jc" Zm Jc 2m 



^ _i [ dx f dy' , 

+s I / Tf{y,x]x] 

JC" 27ii Jc 2m 



Z[J] = 0, (7i 



where the integral contours C, C and C" successively wrap around the negative real axis 
and the singularities in the left half plane. 

In order to derive the VF3 constraints, it is important to examine the explicit forms of 
these integrated operators. First, by using the formulas in the Appendix B, we have 

dy 



Jc2m "^^'^'6Ji 



{z,x,yuxi 

5 



,,,, - + (7.10) 

dJi[y',x,yi,xi,- ■ ■) „=o 

where the second term stands for the polynomial part with respect to y' of the large y' 
expansion of Jy ^ — y The meaning of eq. (|7.1CI| ) is as follows. In terms of the domain 
length /' conjugate to y', the l.h.s. represents the operator with the domain length /' + e, 
where e comes from the procedure of changing the contour /c ^ ' ' ' = l™£^+o P-^oo ^^^^ ' ' 
When /' 7^ 0, the limit e —>■ +0 is smooth, and thus it coincides with the first term of the 
r.h.s. However, for the singular terms supported only at /' = 0, if any, the l.h.s. gives no 
contribution, since lim^^+Q S^'^\e) = (n = 0, 1, 2, ■ ■ ■) in the prescription of the Appendix 
B. The second term is needed to subtract such a contribution. As is seen from the expansions 
( |6.15| ) ^( |6.18| ), for example, c^'s are all zero for I = A,B, 1, and = of, Cg = a^x + a'^xi + 
a'^yi, the others vanish for 1 = 2. Also, for composite operators, the above formula can be 
used. For instance, 



/ -t^z ViV — = = = = = h> y c„ (x). 7.11 

Jc2m '5Ja{z) 5 J,{z,x) 5JA{y')U,{y',x) ^ ' ^ ' 



From the short-domain expansions, c^^(x)'s turn out to be polynomials of x and to vanish 
for n > 4. Thus, we have 

/ — / ^D^{y.y')^ ^ = __^ / ^ J , (7.12) 

Jc" 2m Jc 2m 5Ja{z) 5Ji{z, x) 5JA{y') ■'C" 2vrz 5Ji{y', x) 

and similarly 

/ I Dix x') ^ ^ = I ^ ^ (7 13) 

Jc'2mJc2m ' ' 5 .1 b{z) 5 .h{y , z) Jc" 2m SJeix') 6Ji{y,x')' ^' ' 
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After using the above formulas and doing the integral in the r.h.s. of ( [7. 131) by substituting 
the expanded forms (|6.16| ), ( |6.17] ), the integrated operators are written as 



_i /• dx f dy' , 

s TT- TT-Ti {y;y ,x 

Jc" 2m Jc 2m 



-J— Oi(i)(y)-c02(i)(y) 
dJA\y) oJ dJ 

6 



+s ^(a - ail)— ^— 5i/3 



too ^7/' _ X 

^.JA{y')dy'D-Ay,y')0,{-^){z) 
ioo Zm oJ 



+ (terms containing J/ (/ 7^ A)), 



(7.14) 



-1 f dx P doc I 



s-^ — —Tf{y,x;x') 
JC" 2'Ki Jc 2Txi 



-s ^a^i^a\ji^y'^ — s ^(af/aa'-^r/a + a:^2/3'^-i/3)^?/^ 



^^(«f/3^-7/3 + 4/3)^2/3^ 



2 



-1/ B ///I |1ND -IB II rp 

-s {a^i^a__^i^ + a5/3ji^5/3 - s a_2i^a_ii^i 
-cil + c~\s-'af/,a%)0,i^)iy) 

I B {iv)l 1 \ D 

+s [a - a4/3«-7/3 - ai/3)-y— -5i/3 

+ (terms containing Jj {I A)). (7-15) 

The disk parts of ( 7.14 ) and ( |7.15|) can be interpreted to represent the continuum versions 
of the Schwinger-Dyson equations ( |C.2| ) and ( p.3|) , respectively. In (|7.15|) , thanks to the 
short-domain expansions, the integral J^j ^ijpj^jj'^^^ whose ^/-dependence was difficult 
to see in this form turns out to be essentially a polynomial of y. This property is crucial for 
the derivation of the constraints. 

Also, from the results ( p.l2| ) and ( p.l3| ) in the Appendix D, we can assume the following 
identity holds: 

_^ r dx' r dy' _^ f dx 



c" 2m Jc 2m Jc 2m 6J2{y, x, y', x') 
_i f dx' r dy' _^ f dx 5 



c" 2'Ki Jc 2m Jc 2m 6J2{y', x, y, x') 
_^ r dx' f dx _^ f dy' 6 



c" 2m Jc 2m Jc 2m 6J2{y, x, y', x') 

s-^l ^[ / . (7.16) 

Jc" 2m Jc 2m Jc 2m SJ2iy, x', y', x) 
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In the Appendix D, we have confirmed these equations for some simple cases. This means 
that the results of the spin-flip processes are independent of their orderings. By using these 
identities for the explicit forms of T{^, T^, we can see the symmetry properties: 

dx f dy'^ 



'1^1 ^T,\y-y\x)Z[J] 
Jc" 2m Jc 



dx f dy' , 



c" 2Txi Jc 2ni 



Ji=0 (If^A) 

Tf{y';y,x)Z[J] 



J 1=0 (/^A) 



^1 



Jj=0 (I^A) 



s-'l ^1 ^Tr{y,x';x)Z[J] 
Jc" 2m Jc 2m 



doc f dx 

'l'{^[y,x;x)Z[J\ 

Setting J/ = (/ 7^ A), eqs. ( [7.8|) and (|7.9| ) become respectively 



.-J— + 0',- cO.i^My) + s-\d - a'l,)-^B 
dJAW oJ 6J oJAiy) 



^.JA{y')dy'D.{y,y')0,{-^){z)]Z[J] 



1/3 

0, 



J 1=0 (I^A) 



(7.17) 



(7.18) 



.-l^B A 



+ «2/3)^2/3l/ - S"^(af/3a'"^/3 + 4/3)^5/3 - s'^ a^y^a'^y^T"^ 



s l,a.4/3a_7/3 



c{l + c~h~'a^/,a%)0,ijj){y) 



+ s ^(a - af/3a^7/3 - cti/s) 



SJAiy) 



Bi/3]Z[J] 



Ji=0 (I^A) 



We used ( [7. 141 ), ( [7.15| ) and ( |7.17| ), and O'q is the ^/-independent operator 

On, = -s I / T, (y,x:x). 

° 2 Jc" 2m Jc 2m ^ ' ' ' 



Also, from the once- integrated version of ( |7.1| ), we have 



cOi(-L)(y) + O'' 

dy' J , , 



0, 



Ji=0 (I^A) 



(7.19) 



(7.20) 
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where Oq is a ?/- independent operator introduced as an integration constant. 

The three eqs. ( [7.18|) , ( [7.19|) and (|7.2(]| ) lead to a closed equation of the loop operator 



6'' 



6 



dy' J 



SJA{y) 



+& I i^-JA{y')dy>D^{y,y') 

J —ioo ZTTl 



6 



6 



+ 



52 



+^^st 



too 
ioo 



^5Ja{,z) 5JA{y) 5Ja{,z) 



dy' 



'-too 2Tli 
X Z[J]\,jj=Q (^i^A) 



dy" ~ 



77^ / ^JA{yVA{y")dy,dy.D,{y,y")D^{z,y' 



5 



5Ja{w) 



0. 



The coefficients a's are defined by 



a4 



CS a^l'iO'-'jl'i 

1 + c-is-iaf,,a(")' 



a2 



ai 



4/3"-7/3 

CS Haf/ga" + 02/3) 



1 + c-is-iaf/,a^"^' 



4/3"'-7/3 



CS "[aZ/s^-i/s ~^ ^-2/3^-1/3) 
1 + c~^s~^af^^a[_y^ 

CS ^{a^/sd'Hys + al/s) 
1 + c~^s~^af/^a[_y^ 



ao 



CS ^a^2/3^-l/3 



aA 



— 1/~ B 1 \ 

CS [a — a4/3a_7y3 — a^^^gj 



CS {a — a 



Further, by rescaling as 
6 



SJAiy) 



16 J 5JA{yr 



My) 



9st 



-04 \ 



1/3 



16 ; ^st' 



T 



-a4\-i/3 ~ 



16 

04 



T, 



02 



and putting 



eq. ( p'.21| ) takes the form 

[-IQy^ + lQTy^-yOA-O^ 



a4 



-a4\2/3 



(7.21) 



00, 



dy' , 



2 g|2 

'5JA{y) 



/too f^y' _ 



5 



5 



+ 



52 



^6Ja{z) 6JA{y) SJa{z) 
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+^st / [ 7r-JAiy')JAiy")dy>dy>>D^{y,y")D^{z,y" 

J —too ZTll J— ^oo Z7C1 



-ioo 2Tii J-ioo 2Tci 5Ja{w) 
X Z[J] |j,=o (W) = 0. (7.22) 

Let us now confirm that this equation represents the W3 constraints. In order to do so, 
we have to identify the disk and cyhnder singular parts. For the disk amphtude w{y), the 
Schwinger-Dyson equation obtained from ( [7.22| ) is 

w{y)^ - {Oo)ow{y) - IQy^ + 16 V - 2/((^a)o - (0i)o = 0, (7.23) 

where we consider the solution regular except on the negative real axis. Note that this 
condition does not determine the solution uniquely. In fact, we find two solutions: 

i) w{y) = 2'/'{y-^){y + 3^y/^ 

(Co)o = 0, (Oa)o = 128 j , {O^)o = --T\ (7.24) 

ii) w{y) = {y + Vy^Y'' + (2/ - Vy^Y'\ 

{Oo)o = 3T^/^ (Oa)o = 0, (Oi)o = 2T2. (7.25) 

It is the solution ii) that reproduces the matrix model result. By repeating the argument 
in the Appendix C for deriving ( p.l4| ), without using the Z2 symmetry A ^ B, we can 
see that (Ca)o is proportional to the next leading order (0(a^)) of the universal part of 
the quantity -^{tr^A — -B))o- This implies that the solution i) spontaneously breaks the Z2 
symmetry. Here we only consider the Z2 symmetric solution ii) and leave the case i) as a 
future problem. 

For the cylinder amplitude w{y,yi), the Schwinger-Dyson equation is derived from the 
lowest order of g^j^ in the Ja(z/i) derivative of ( [7.22| ): 

-y{OA^A{yi))o - {Oi^A{yi))o - {Oo^A{yi))ow(y) 
+ {3w{yf~{Oo)o)w{y,yi) 

+ 9lt9yi:i^^iMy)'-w{y)w{y,)-w{yif) = 0. (7.26) 
If we take the solution ii) as the disk amplitude, note that it satisfies 

3w{yf - {Oo)o = at y = 0, ±W ^. (7.27) 



The disk amplitude is regular except the region y < —yT on the real axis. By assuming 
the same property for the cylinder amplitude, eq. (|7.26|) can be solved, by using the similar 
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argument as for the case ( |B.7|) . The result is 

4 1 



^Stn 



X 



9 f{y, yi)9{y, yi) {y + Vy^y/^ + {y - Vy^^f' + 

{yi + ^yl-Tf'' 

{yi + \lyl - TY" + {yi - ^y\ - Tf'^ + 



X 



f{y,yi) 



9{y,yi) 



(7.28) 



where 



+(?/ - ^Jf^y^'iyi - M^Y" + T^/^ 

{y + 77^)^/^ + (yi + ^Ivl^f" 



9iy,yi) 



From ( |7.25| ) and (|7.28|) , the singular parts can be found as 



w 



w 



smg 



smg 



(y) 



iy,yi) 



2ny'"-\y-'% 



alt 



(7.29) 



9 {y - yif 



The connected correlation functions are expanded by the local operator insertions ga^^ 



as 



{^A{yi)) 

(I>a(i/i)$a(|/2)) 

{^A{yi)---^A{yK)) 
9^''\yi,---,yn) 



w'^Hy,) + 9^'\yi), 

w'''''S{y,,y2)+9^'\yuy2), 

9^''\yi,---,yK) {K>3), 



(7.31) 



where a^'s run over the positive integers +1/3 and +2/3, i.e. 1/3, 2/3, 4/3, 5/3, 7/3, ■ ■ ■ . 
Using ( |7.29| )^( |7.31| ), we expand ( [7. 221 ) similarly as in the argument of the Virasoro con- 



straints in the Appendix B. From here the analysis is parallel to the reference ||T9[, where 
the constraints were explicitly derived from the two-matrix model for the first time. So, 
we show only the results. For Oq, Oa, Oi insertions. 



(Oo) 



3 ■ 2^/3 



9l/3, 
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"J ai 

((^aI'a(?/i)) = 42«/3|yr^/'+3-2«/3$:52A..yr-\ 

ai 

(Oi) = yT^ + 3 ■ 28/3^5/3, 
(Oil'A(yi)) = glt^'^'^yl^' + 3 • 28/3 J: g,/s,^,yr'-\ 

"J ai 



= 3-28/3 ^ ^^^^^^^^ 

{K > 2). (7.32) 



This means that Oq, Oa, Oi are expressed in terms of the loop operator 

OJA 

dy ^.1/3 5 



^ 3.24/3 / ^^1/3 ^ 

= 3.28/3/ ^^V3^^ 

Jc 2m 5JA{y) 

O, = 1^^2 + 3-28/3 / ^^5/3^ 

3 Jc 2m 



SJAiyY 

as operators acting on Z[J]|jj-=o {i^a)- 

For the other contributions, by introducing the generating function 

00 "1/3 "2/3 

V — ^ i\/3 / 2/3 

In = — — - - - 51/3,. . .,1/3,2/3, . ,2/3,.. , (7.33) 

n„=0 ^1/3' ^2/3! ^ ^ ^ ' 

they can be expressed as 



"1/3 "2/3 



L„Z(/x) = (n>-l), 
Wn'Z{pi) = (n' > -2), (7.34) 



where 
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' ■ ' [ai^, - 13^,) ^'h^ """as: + 27- 



(l > 1), 



a//8/3 3 dlJ,2/A ) 
^2 



3 r;^' ^/^ 



'a+a'-2 



+4^/^4/3(/il/3 - ^st 3 • + ^^t|/^2/3, 



= 3-28/M- -T- — 

\ rill-, 1 X riiir- /r. / 



1^9//ii/3 3"^ 9/^5/3 J ^ ^ dii\ 

+5'st E E ^ — ^ — 

+.s\^(/^i/3-%^3.2^/3T)3, 



■2/3 



+3 • 2^/=^ ( E E I 

+ E 5' 

9si E (()/Xa+m+7/3 3 (?/fa+m+l/3 ^ 

52 

+^St E E ^ ^ 

/3+/3'=a+m ^^I^P^I^P' 
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1 d 
+5'st o aa'fiaf^a'T. (m > 0). (7.36) 



This is nothing but the W^s constraints. 



8. Conclusion 

Let us first summarize what we have done. We have started our paper by discussing the na- 
ture of PSFTs from the view point of the stochastic quantization of the matrix models. Then, 
we have presented detailed derivations of the stochastic Hamiltonians in the double- scaling 
limit from the matrix model, and investigated the infinite algebras of the Schwinger-Dyson 
operators appearing in the Hamiltonians. We have also checked that the algebras contain the 
Virasoro (one-matrix model) and W^, algebras (two- matrix model), as they should. Proofs 
of some of the crucial formulas have not been completed, because of technical complexity. It 
is therefore desirable to develop more powerful methods of treating the double-scaling limit 
for general target spaces. 

After these calculations, we have to reconsider the questions raised in the earlier sections 
of the present paper. Perhaps, one of the most important lessons of our work is that the 
structure of the general splitting and merging interactions of string fields with arbitrary 
matter configuration is not affected by the mixing of the string-field components in taking 
the scaling limit which is defined for a specific background. This seems to imply that the 
structure of these terms is completely independent of the backgrounds. Recalling the general 
discussion in section 2, we realize that the purely cubic Hamiltonian of the matrix model with 
most general source terms and no bare action already captures the structure of the continuum 



Hamiltonian in a background-independent way. This conforms to earlier suggestions |£l 
and points to an intriguing possibility of formulating a background-independent string field 
theory, encompassing critical strings, by starting from general matrix integrals with infinite 
number of different matrices. For the case c < 1, a related idea has already been discussed 



in ref. [15 



Before pursuing such possibilities, there remains, however, many important issues to 
be solved. Asides the problems mentioned in the Introduction, what is needed to make 
further progress is a deeper understanding on how to extract real space-time picture of the 
string theory from matrix models, since matrix models apparently miss some important 
characteristics p2|, E3| of the string dynamics. 
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Appendix 



A k = 3 critical point and the disk amplitude 

Let us begin from deriving the non-even critical potential for k = 3. It is sufficient to recall 
the well known formulas in the method of orthogonal polynomials. For the coefficients 5* 
and R in the recursion equation for the orthogonal polynomials Pn{^) = A" + ■ ■ ■, 

XPniX) = P„+l(A) + 5„P„(A) + RnPn-l{X), 

we have a set of equations in terms of the potential V, 

= i^-V'iz + S + -), (A.l) 

in the sphere limit ~ — oo, where x = ^, S = S{x) ~ Sn, R = R{x) ~ Rn- Eqs. (|A.1| ) 
and ( |A.2| ) give the relation that implicitly determines the function R{x), of the form 

X = W{R). 

Since the free energy is determined by R as 

N-l 

In Z ~ ^ (A^ - n) In (A.3) 

n=l 

the following behavior of the W{R), 

W{R) = l~const.{l- Rf (A.4) 
as X —>■ 1, R 1, leads to the k = 3 criticality of the free energy 

lnZ~(l-|r/3. 

This shows that the minimal order of the k = 3 critical potential is four. After lengthy 
calculations, we find that the potential 

V(M) = ^(^M2 + + — M^) (A.5) 

^ ^ A^^ 2 3 20 ^ ^ ' 
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with (|3.19| ) and (|3.20|) satisfies all of the above conditions. 

Next let us derive the disk amplitude in the sphere approximation. By using the method 
of ref. the disk amplitude in the large N limit is given by 

(<^>(C))o = ^v^'(C) + Ij^i-l^ + + i?)\/(C-M(C-&-), (A.6) 



where 



A = -(73-^(&+ + fe-), 



B = -92-^ib++b-)-^{b^+b^r-^{b^-b^r. 

By introducing the variable 

z = 10(73 + 3(6+ + (A.7) 
the end points of the cut b_ < ( < b+ are determined by the equations: 

- + 36z^ - 256z - 720 + 81920^"^ = 3888^, (A.8) 
4 (3 

I. 2 2 80 2560 1 ^, 

The disk amplitude with a microscopic loop can be read off from the coefficient of in 
the large ( expansion of eq. (|A.6|) , 

(^trM)^ = (&+ - b-fMb^ + 6-) + gs{b^ + b^f 

+1 - b^f + ^(6+ + b.f + + _ (A.IO) 
In the scaling limit 

^ = l-aH (A.ll) 
(a means a lattice spacing), eq. ( |A.8|) is iteratively solved as 

1 + a-t^l^ + a^—t^l'^ - a^^^t 
4 64 1024 

Then h+ and 6_ are 



o 49 

3 24 384 

165888 ^ ' ^ ' 

I^^i^ + a^lt + a^iit^/s + o(a^). (A.14) 
3 16 256 ^ ^ ^ ' 
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Also, it can be seen that ( should be tuned to the critical value at b-^ 

C = C*(l + a?/), C = ^ . 



Substituting ( [7mD ^ ([X:i5|) into and (pTTOD , we have 



where 



($(0)0 = lv'iC) + a'/My) + Oia'/'), 

/I ,A 32 + 25^3 o4 ,3 

(— trM) = — - a^-t + a^-t^^^ + 0(g 

\N /o 15 5 4 ^ 



(A.15) 

(A.16) 
(A.17) 

(A.18) 



with T = (2C')=^t 



B Expression for the local operator Oq 



We give the proof of eqs. ( [4 .121) and (|4.13|) using the Schwinger-Dyson equations. We first 
rewrite the once-integrated Schwinger-Dyson equation 



Tiy)Z[J] = 
into the relations among the local operators. Here, 

T{y)=To{y)+ f dy' piy') + O,, 
Jo 



(B.l) 



(B.2) 



the y- independent operator Oi appears as an integration constant. 

Now to correctly extract the local operators, we must identify and subtract the singular 
parts in the correlation functions 



6^ \nZ[J] 



dJ{yi) ■ --SJiyK) 



{^y^)■■■^yK)). 



(B.3) 



j=o 



It appears only in the disk and cylinder amplitudes. For the disk, from the large y expansion 

of 



we see that the first two terms of 



correspond to the singular part 



^sing(^) = y'/' + i-Ty-"\ 
lb 



(B.5) 
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because it is known that the disk one-point function of the local operator behaves as T'^l'^+^ 
(A > 0) from the analysis of the continuum theory |]TB|. 

For the cylinder amplitude, we start with a derivation of the amplitude w{y, yi) from the 
lowest order of g^i in the Schwinger-Dyson equation 



^ -ny)Z[J] 



that is, 

2wiy)wiy,yi) + gl^dyMy,y,)wiz) - y(Oo$(yi))o - (0i<i>(yi))o = 0. (B.7) 

It can be easily solved by noting the fact that w{y) given in ( |3.37| ) has single zeroes at 
y = a,a, where a, a are the solutions of the quadratic equation 

- \T"'y + It'" = 0, 



explicitly 



a 



- (B.8) 

a 4 



By setting y = a or a in ( p.7| ), the first term vanishes, fj\ then ((9o$(l/i))o) {Oi^{yi))o are 
determined as 



(Oo$(yi))o = 49,,Vyi + TV3, (B.9) 



{OMyi))o = gltdyAyi'l;T'")\/yi + T'"- (B.io) 



Substituting these into ( |B.7| ), we obtain 



2 

Hy, yi) = ^ , \ , ^ , (B.ii) 



Since the cylinder amplitude of the local operators behaves as 7"1/3+Ai+A2 (^^^^ A2 > 0), the 
singular part is identified to be 

4 - yiY \\l yi \l y J 

Thus, the connected correlation functions are expressed as 

i^y,)) = nJ'''''Hy^)+g^'\y,), 

{Hyimy2)) = w''''Hy„y2)+9^'\yi,y2), 

{^y^)■■■^yK)) = g^^'Kyi, ■ ■ ■ ,yK) iK>3), (B.13) 



Here, we assumed that w{y,yi) does not have any poles aX y — a, a. This is justified, since is 
regular except in the negative real axis as is seen from its definition. 
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where g^"'\yi, ■ ■ ■ , Vn) is the part interpreted as local operator insertions, and it is expanded 
by the correlators among local operators gai,-,a„ 



(B.14) 



ai, •••,«„ 



where a^'s run over the positive half odd integers 1/2, 3/2, 5/2, • • •. Using ( [B.5| ), ( |B.12D ^( |B.14| ), 
we expand the Schwinger-Dyson equations 



6 



K 



5J{yi) ■ --dJiyK) 



ny)Z{J] 



(i^ = 0,1,2, 



(B.15) 



j=o 



and perform the similar analysis as in ref. [l^l. For example, the result for = is 



-y{0,)-{0,)+[-Tj y-' + gitY^y- 



^ lo 



+ Y^ig^g^' +ga,a')y 



-a-a'-2 



Here, we use the Greek indices a, a', (3,13' for the positive half odd integers. From the first 
two powers of the large y, {Oq), {Oi) are determined as 

(C»o) = 2^?i/2, (Oi) = 2^73/2. 

Performing similar analysis for i^' > 1 in (p.l5| ), we obtain for Oq, Oi insertions 

(Co) = 2^71/2, 



-.2 -"^^-l/S 

ai 

{0,^y^)---^{yK)) = 2 ^?iAa„...,a,.l/r^-'---Z/x"""\ 



(Oo$(2/i)) = glt7^yi+^j:9i/2,a,yi'''"\ 



(Ol) = 2(?3/2, 

{OMVi)) = gltlyl^' + 2Y^gy,^^,y^-^~\ 

{pMyi)---~^{VK)) = 2 9V2,a,,-,a^yl'''~'---yK''"-' {K>2). (B.16) 

This shows that {Of)^{yi)) and have the singular parts gl^^yi^^'^ and 

respectively. 

The other powers of y give the Virasoro constraints. For K = 0, we have 

5 

2(^i+7/2 + T^^^^+i/2) + J2 idpgts' + gp,p') 

I3+I3'=l 



5_ ^2 
16' 



+9ltY^i,o+{-^T) 6,., = (/ = -l,0,l,---), 
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where g with negative indices is understood as zero. For general K, by introducing the 
generating function 

^ "1/2 "3/2 
X — ^ '^1/2 ''3/2 

\iaZ{fi)= 2^ r ••■^1/2,... ,1/2,3/2,.. ,3/2,..., (B.17) 

nt^O ^1/2 ! ^3/2 ! ^ . ' ^ . ' 

"1/2 "3/2 

the relations can be expressed in the form 

L,Z(/i) = (ri>-l), 

= '5£;^+^^tE«/^«^ + (^^+4i/^l/^)^ 

9 5 (9 (9 1 

(/>1). (B.18) 



Note the operation of Oq (Oi) on 2'(yu) is expressed as the local operator insertion 2- 



d 



'■9^1/2 

d \ 



Next, we notice that as an operator acting on Z{^), is expanded by the local operator 
d . 

s d 
— — = - — -u^"^^ as acting on Ziu). 

Furthermore, since the partition function Z[J] is related to the Z{^) through the rescaling 
dy ?....,sing. ^ _^ 1 [dyi f ^T(,,^T(,,^,M 



z[J]=exp[j ^^Ay)w'''''Hy) + ^J ^Jiyi)Jiy2)w^^^^ (b.i9) 

we see that 

-i- = w'^^Hy) + f p:J{y,)w'^Hy, y,) + ^ -^y—^ (B.20) 
6J{y) J 2m ^ o^ia 

as acting on Z[J]. 

After these preparations, we can now prove the eqs. ([4.1 2|) and (|4.13| ). First, we consider 
( [4.12|) . The integral along the contour C is defined by the analytic continuation using the 
Beta function. For example, 

dy A dy . .^.^ , dy 



Jc 2m J-oD 2m Jo 2m 

' , A sinvrA , . 

/ dyy^ = B{A + 1,-A-1) 

Jo TT 



TT JO 
1 1 



A + ir(o)' 
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For A —1 this is zero, and for A = —1, because the integrand has no cut, the contour can 
be deformed as encirchng the origin 

Jc 2m Jo 2m 

We summarize these into 

/ = Sa,-i. (B.21) 

Jc 2m 

From this result, in the case that the pole y = —yi exists inside the contour, we obtain 

Jc 2m y + yi ^ Jc 2m 



\ i-yi^ A = 0,1,2, 
I otherwise. 

Also, when the pole y = yi is outside the contour, 

f dy y^ _ [ A = 0,l,2,--- 

Jc 2m y -yi \ -yt otherwise. 

By using these formulas (|B.21| ) and (|B.23|) , we can see that 



;B.22) 



(B.23) 



Jc 2m ^ Ojla OjJLii 



/2 

y'/V'^Hy) = = {Oof''^ 



Ic 2Tii 

2 ^y'/'w^^'^Hmyr) = gltlyl'^' = {OoHy^)f^^, 

which just mean that Oo is written as in (|4.12|) . 

Next, we verify the second equahty ( ^.13| ). For this purpose, it is sufficient to show that 

lim r ^e^-vy^ = I ^y^ (B.24) 
£^+oj-joo 27ri Jc 27ri 

for arbitrary A, because the formulas such as ( p.22|) and ( B.23| ) can be derived from ( p.21|) . 
We consider the following integral in the case A ^ Z 



iOO 



—y^ = -cos / rfxx^ = -cos B(A + 1,-A-1) 

2m TT 2 JO vr 2 



too 

1 1 



r(0) A + 12sin 



nA 



0. 
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This implies that for A ^ Z, 

dy ev A ^ ^" r°° dy 



lim 



J-ioo 2-Ki £^+0,^ n! J -too 2ni 



Also, when A = — 1,— 2,— 3,---, the contour can be deformed to the circle enclosing the 
origin, 

lim r ^e'Vy^ = hm / ^e^^y^ 

s\A\-l 

Further, for A = 0, 1, 2, ■ ■ -, the derivative of the 5-function aX y = e appears. In our 
prescription, since S^'^^e) = for finite e, the limit e +0 is also zero 

dy 



lim n ^e^yy^ = lim 6^^\e) = 0. 

:^+0 J-ioo 2m 



These results are summarized into 

lim r ^e^yy^ = Sa,-i. (B.26) 

Comparing this with ( p.21]) , we conclude that eq. ([4.13| ) holds. Similarly, it is easy to see 
that the operator Oi is expressed as 

Jc27n^ SJ(yy ^ ' 

= 2 lim r ^e^yy^/^^—. (B.28) 

e^+Qj-too2m 5J{y) 



C Disk Amplitudes in the Two-Matrix Model 

Here, we obtain various disk amplitudes (genus zero one-point functions) in the two-matrix 
model by using the continuum limit of the Schwinger-Dyson equations which give the re- 
lations among them. Some of the disk amplitudes before taking the continuum limit have 
been obtained by Staudacher . We will extend his results considerably and give detailed 
forms of the continuum disk amplitudes that have not appeared in the literatures. 

We introduce the following notations for the disk amplitudes (some of which are borrowed 
from 13): 
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^(0 = /^tr ^ 



W{a) = /Itr 







c 


-A 




1 


a 


- B 




1 


c 


-A 




1 


c 


-A 




1 















H'<-'(C..:..-.a...) ^ (Itr^^...^^)^ (.^1.2....). 

a J w(c) 

Let us first start from the disk ampfitude with the simplest spin configuration on the loop. 
(The spins on the loop are all ^4-.) It is obtained by combining the following three Schwinger- 
Dyson equations: 

(C-K'WC) = cW,{0 + w(0' + i-9(C + w,), (c.i) 

iC-gOWiiO = cW2iO + WiOW^{0 + W^-g{CW^ + wi'}), (C.2) 
W^iO-gW^iO = cCWiC)-c. (C.3) 

We can eliminate H^i(C) and W2{C), and have a cubic equation of W{(): 

W{Cf + aiW{Cf + a2W{C) + as^0, (C.4) 



a2 = (C-Ky--(C-K') + (--^)C + l-^W^i, 
as = (-l + (7iyi + (?C)(C-K') + (l-3c + cK)W^i 



yVF3-^+-(l-c2)-cC, 



(2) 

where in order to eliminate W2 and { we used 



Wi - gW2 = cWi, W2 - gWs = cW^/? + 1. 
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The expressions of Wi and Ws are evaluated by the orthogonal polynomial method ||T^ 
as follows: 

Wi = ^[3p^-6cp3-2(l-2c)p2-2c(l-2c)V^ 

+32/- (1 -2c + 4c^)(l -2c)], (C.5) 

w, = Y^:^[-i6(/-i) + 90c(/-i) 

+ (80(1 - 2c) - ^c^) (p^ - 1) 

+ (-64 - 94c + 380c2 + 60c3)(p^ - 1) 

+ (-48 + 336c - 333c^ - 150c^ - 54c^)(p2 - 1) 

+2(1 - 2c) (32 - 41c - 14c=^ - 66c^)(p - 1) 

-6c(l - 2c)^(7 - 14c - 2c^){p-^ - 1) 

+c(l - 2c)2(16 - 21c - 6c^ - 6c^){p-^ - 1) 

-4c3(l - 2c)3(p-3 - 1) - ^c^l - 2c)\p-' - 1)], (C.6) 
where p is implicitly determined by 

/ = -^[4p3-9cp^-4(l-2c)p + 2c(l-2c + 2c2) 

-c(l-2c)V1- (C.7) 
In the continuum limit, expanding g and p about the critical points 

9* = \/l0c3, P* = 3c, (c, = ^ 27^^ ^' 
eq. ( |U.7| ) can be solved iteratively: 

P = P. + a'/'lp^{5ty/' + a'/'^o^{50' 
3 36 

where g is expanded as g = g^:{l — a^t). 

Substituting this into eqs. (|C.5|) and ( |C.6|) , we have Wi and in the expanded form: 

Wi = W^'"' + Wi, (C.9) 

^ -8pt + 3{2g,y , -136pt + 27{2g,y 
' 3(2(7,)3 27{2g,y 
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27(2^*)3' ' 81(2^*)3 

!7p^ + 97 
972(2^,) 



= iy3"°>^ + #3, (C.IO) 

^ 32(420 - 839p,)p5 , 160(252 - 611p,)p5 
=^ 729(2(7.)^ 729(2(?,)5 

M/ _ 8/3 320p^ rc;f^4/3 I ^10/3 160P* rp;f^5/3 

729(2(7*)^ 

where we denoted the non-universal pieces by Wi°'^, W^°'^ and the universal ones which give 
the continuum limit by Wi, W3. 

Now, we shall evaluate W{() in the continuum limit. Shifting W{() as 

w{c) = -^ + wiO, (c.ii) 

eq. ( |U.4| ) becomes 

W{Cf-lA2W{O-^A, = (C.12) 

where 

Ai = 9aia2 — 2a^ — 27a3, A2 = a1 — 802. 
Then the critical point of ( denoted by is determined by 

A,l = A2U = 0, (C.13) 

where |* means that g, Wi and W3 are set to the critical values. It turns out that eq. ( |C.13| ) 
gives a cubic equation of P,,, and its solution is threefold: P^, = ^4^- 



After substituting ^ = P^,(l + ay) into ( |C.12D and expanding with respect to a, ( |C.12| ) 
becomes 

where c is fixed to be the critical value c*, s is the irrational number s = 2 + -\/7, and the 
rescaled variable T = is introduced. 
The solution of (^^) is 

„l/2 „4/3 ^ ^ 



56 



which gives the universal part of the disk amphtude. 
Also, the non-universal part W'^°'^{() is 

w^^^O = -T = -f + l(C - gC')- (C.16) 

C.2 W^{C),W2{() 

The amplitude Wi{Q (W^2(C)) represents the configuration that the spins on the loop all 
align A- except a small i?-domain consisting of a single spin (two spins). 
From (ICH ), 

w,{C) = ^^[{c-gc^)w^°-{c)-w^°^{0'-i + g{c + wr)] 

+l[C-ge-2W-"-iO]W{0 

+^[-WiCy + 9W,]. (C.17) 

We identify the universal and non- universal parts as follows. If there are polynomials of y and 
T, they are non-universal. Also, if there are amplitudes, with the spin configurations simpler 
than that of Wi{Q, multiplied by polynomials of y and T, they are non- universal. After 
these identifications, the remaining terms are universal. By using this rule, the universal 
and non-universal parts, denoted by W^i(C) and Wi°'^{() respectively, are determined as 

w,ic) = wrio + mic), (C.18) 
+}-[(- gc'-2w^"-{omo, (C.19) 



a 



8/3 _ 



S 



8/3 



40 ■ 22/3 I 



= a"'^-^,wM + 0{a'). (C.20) 
By a similar manipulation for eq. (|C.2 ) using 

we have 



8/3 



1 C 1 , Ox 
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+ 



w{0 



W2{0 



9 

--W{C)W,{C) + 0{a'''') 
c 

r.4 



a 



leovioc 

160a/10c 



(16/ - IQTy' + 2T') + 0{a'"'') 



C.3 W^^\C,cr),W^'\Ci,criX2,cr2) 

To discuss the higher disk amphtudes, the recursion equation for W'^'^^^ given by Staudacher 
(eq. (20) in ref. (H) is useful: 

X {Dq(Qi, g,) Q[g - - w^(gi) - w{Qk)] - 1 + ^^Pi) 



+- EPq(Qi' Q/)w^^''''^(^2, Q2, ■ ■ ■ , P/, Q)] 



fc-1 



1=2 

X [/^Q(g^ gfc)w^(''-''H^/+i, Q^+i, • ■ ■ , Pk, Q)] 

+cDpiP,, P2)W^''-^\P, g2, ■ ■ ■ , Pk, Qk) 

-j2W^''~'\PuQi,---,Pi-i,Qi-i) 

1=2 

X Dp{P„ P,)l^(2^+2-20(p^ g^^ . . . ^ p^^ g^)|. 

Let us consider the continuum hmit for the cases k = 1 and k = 2: 

{l-gOW{a)-cW{C)-gW,{a) 



iy(^)(Ci,ai,C2,cx2) = 



C-gC-ca- W{0 
1 



(C.21) 
(C.22) 



C2 - ^?CI - ca, - WiC2) 



58 



X{(C - W^'\Cl. ^2))/^c(Cl, C2)W^'\C. ^l) 

c 

+ -w{Ci)}. 



(C.23) 



For k = 1, putting ( = P^,(l + ay), a = P^,(l + ax) and expanding with respect to a, 
( ^2%) becomes 



p^(2)non(^^^) = c{l-as{y + x)) + VlO^{W{() + W{a)), 



(C.24) 
(C.25) 



+a 



,-l^(C)H^(cr)-H^(C)' + cW^i(a)) 

^1/2^2 



(y + x) 



cs 



[160(/ + x^) + 80{yx^ + y^x) - 40?/V 



120 {y + x)' 



+ 12{s - lO)T{y^ + x^) + 24(s - 5)Tyx - 120T{y^ + x^ + yx) 
+ l2sT{y + x)^ + 15T2] 
+0{a'/'') 



5/3 



4 • 22/3 



where 



w^'^\y, x) 



-w{y)'^ - w{y)w{x) - ^(a;)' + 3T^/^ 
y + x ' 



(C.26) 



(C.27) 



and we calculated W^'^\C, a) up to the next leading order since it will be necessary to obtain 
W^''\Cu(^iX2,cT2) below. 

Repeating the same procedure for k = 2, we have 



W^'\CU CTU C2, ^2) = ^^('^'^""(Cl, ^1, C2, 0^2) + W^'\Cl, CTl, C2, CTs) 



(C.28) 



W(^)'^°^(Ci,ai,C2,a2) = 10c2-10c[Dc(Ci,C2)l^(C) + i^.(^i,^2)l^(a)] 

-VWc[D^{Cu C2){w^'\C, <yi) + w^'\C. 0-2)) 

+D<,(ai, a2)(iy(')(Ci, ^) + W^'\C2, a))], (C.29) 



l^^'HCi, ^1, C2, ^2) = a^^W(yi, xi, 1/2, X2) + 0(a^/3) 
Sc's 1 



(z/i - y2)(a;i - X2) 



-^(yi - l/2)(a;i - X2){yi + y2 + xi+ X2) 
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+ ^{w{yi) + w{x2) + 2w{y2) + 2w{xi))w'^^\yi, X2) 

+^(^(2/2) + w{xi) + 2w{yi) + 2w{x2))w^'^\y2, xi) 

-^(^(1/2) + w{x2) + 2w{yi) + 2w{xi))w^'^\y2, X2) 
+0{a^'''). (C.30) 

We note that the above results have the symmetry under the cychc permutation of variables 

Ci -> 0-2, 0-2^(2, C2^cri, CTi^Ci, 

which should be satisfied by the definition of the amplitude. In general, however, in the 
expressions such as eq. ( C.23 ) the symmetry is totally obscure. Obtaining amplitudes 



with correct symmetry property thus constitutes a quite nontrivial consistency check for the 
continuum results. 

C.4 vri(Ci; (2,0^2) 

From the following Schwinger-Dyson equation 



= / d'^'Ad'^'B E ^ ftr f -^t--^^ 



we have 

i^i(Ci;C2,cr2) = -\[Dc{CiX2){C-ge-w{Ci)-w{C2))w^^\c,<y2) 

+ 9W{a2)]. (C.31) 
By a similar calculation as before, we find 

w^i(Ci; C2, CJ2) = wriCi; (2, 0^2) + C2, (72), (C.32) 



w-^'^iC; C2, CT2) = -D^iCu C2)(C - gOil - 7r(C + ^2 - 2P.)) 

-I * 

+^ g - ^(0 - .CD - ^(C2 - gc!)) + . - - 

-^{W{C,) + W{C2))--Wia2) 
P* c 
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-VWcD^iCu C2)w^i(C) - -D^iCi, C2)((C - 9eW^'\C, ^2)) 

~~c ^^'^ ^^^'0 ^2)^^'nC, ^2), (C.33) 



W^i(Ci;C2,cr2) 



a^j^{w{yi) + w{y2))Dy{yi, y2)w'-^\y, X2) 
+0{a'/''). 



(C.34) 



These results show the following scaling behavior for the general disk amplitudes 
iy(2fc)(Ci,ai,---,a,afc): 

W^''\Ci, (Ti, ■ ■ ■ , (k, (Tk) = J-"^V^\y^, xi, ■ ■ ■ , yk, x^), (C.35) 

which is consistent with the analysis of the boundary conformal field theory 0. An argument 
for this is as follows: The gravitationally dressed spin operator exists at the boundary of 
domains and its dimension is [y]^'''^. This is derived by considering the gravitational dressing 
of the spin operator, whose dimension is [y]^''^, in the boundary conformal field theory in fiat 
space [^. In eq. ( P.35|) , increasing k by one unit corresponds to adding the two domains. 
Clearly, the boundaries of domains are also increased by two. Then, the dimension of w*^^'^^ 
is changed by a factor 

where 2 ■ | comes from the dressed spin operators at the two boundaries, and 2 ■ (—1) from 
the two domains. This coincides with ( |(J.35| ). 



D Continuum Spin-Flip Operator 

In the matrix model before taking the scaling limit, a domain consisting of only a single 
fiipped spin can be obtained as an integral of a general domain, 

A^VC-^ / J 27ii N \C~ Aa- B J ^ ' 

Let us construct the continuum version of this operation. We can do this by deriving the 
relation between the universal parts of the both sides in ( |D.lj ). 

First, let us consider W^'^\C,,a) and Wi{C,). Comparing eqs. ( p.20|) and ( p.27| ), we have 

/ ^.w^'\y,x) = -w{yf + T^'^ = w,{y) - 2^^/^ (D.2) 
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where the contour C surrounds the negative real axis and the pole x = —y. The calculation 
can be performed by using the formula ( [B.22| ) after expanding the numerator of w^'^\y,x) 
with respect to the large x. In such a calculation, we assume that the unintegrated variable 
y is outside the contour, and —y is inside. By including the overall factors, it is rewritten as 

i^aW('\C, a) = W,{0 - a«/3^^r^/3 + 0{a') (D.3) 
J 2m 40 

where the integral symbol §-^cr is used in the sense of 

/"da p2 Z' 
27ri * Jc 2711 

Next, for W'''^\C,i, ai, ^2, (^2) and Wi{C,i] (2, (^2), we use the formulas: 

dXl 1 X? Xn 



(a ^ Z), (D.4) 

c 2ni xi - X2X1 + yi yi + X2 

£L = („ = 0.1.2....), (D,5) 

Jc 2m xi — X2X1 + yi yi + X2 

which are derived from the formulas in the Appendix B, where we regard again that the 
unintegrated variables X2,yi are outside the contour, and —X2,—yi are inside. After some 
calculations, we have 

j^—.w^^\yi,xi,y2,X2) = {w{yi) +w{y2))Dy{yi,y2)w^'^\y,X2), (D.6) 

s-^/^ai#(^)(Ci, ai, C2, ^2) = #i(Ci; C2, ^2) + 0{a'l'). (D.7) 

We can use this method also for a domain consisting of two flipped spins. For a prepa- 
ration, we shall compute the amplitudes 1^2(0; C2, C2) and Vri(Ci; C2, 0'2, Cs; "^s)- From the 
analysis of the Schwinger-Dyson equations similar in the Appendix C, we obtain 

W^2(Ci; C2, (J2) = a^^'^ 320.21/3 ^2^' 

W2{.yi]y2,x2) 

= {-w{y,f~w{y,)w{y2)-w{y2f + ?,T^''')Dy{y,,y2)w^^\y,X2), (D.8) 

5l^2_g4/3 

W^l(Ci;C2,Cr2,C3,Cr3) = ^^'^ . 2I/3 ^1 (^1' ^2, 3^2, Z/S, 2:3) , 

wi{yi]y2,X2,y^,x:i) = iw{yi) + w{y2))Dy{yi,y2)w'^^\y,X2,y3,X3) 

-Dyiy2, y3)w^'\y, X2)Dyiy3, y,)w^'\y, X3). (D.9) 
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Then it is easy to see that the following formulas hold: 

r dx2 I dy2 , _ f dx2 f dy2 

Jci 27ri Jc 2-Ki ' ' JCi 271% Jc 2ni ' ' 

f dx2 f dxi , 

c 2m 02^1""'^^'' ''''''' 
dx2 f dxi 



Ci 2m Jc 2m 



W2 



(l/i)+2T^/3^(yi), (D.IO) 



where the contour Ci wraps around the contour C. Moreover, after a straightforward but 
lengthy calculation, we can show that 



f dX2 f dy2 

dx2 f dy2 



2m Jc 2,,^^(^2'^3,X3,yi,x,) 
f dx2 f dxi . , 
Jci 2m Jc 2m 

dx2 f dxi , , 
wi{yi,X2;xi,ys,xs) 



Ci 2m Jc 2m 

W2{yi;y3,xs)-2T'/^Dy{y,,ys)w^'\y,xs). (D.ll) 



Now by taking the overall factors into account, ( P.10|) and (|D.11| ) are rewritten as 



= W2{(,) + a'/'^-^T'/'WiC) + 0{a''/% (D.12) 
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ol/3„8/3 

= W,{Ci: C3, as) - a'''-^T'I^D,{Ci. C,^)W^'\C, ^3) + 0{a'''% (D.13) 
where the integral symbols are used in the sense of 

/rfcTs 2 f ydC2 „ /■ dy2 Jdai f dxi 

f7r-(^2 = P,a — , f— = P,a —, f— = P,a —. 
J It\i Jci I'm J 2,711 Jc 2,711 J 27i% Jc 27i% 

These formulas show us that the domain consisting of two flipped spins is constructed by 
shrinking the domain between two microscopic domains consisting of a single flipped spin. 
Indeed, by substituting eq. ( |D.7| ), the first line of ( |D.12| ) becomes 

Thus the symbol s~^§^ai is the single-spin flip operator, while shrinks the domain (^2 
to nothing in conformity with the original matrix-model operation. 



E Commutativity of the Mixing Matrix with Sphtting and Merging Pro- 
cesses 



Here, we present the calculations which leads to eq. (|5.28| ) 



6J V 6J, 



M 



,5J 5.1, 



for the first several components I = A, B,l,2, and eq. ( |5.29D 



K,L 



for (/, J) = {A, A), {B, B), {A, 1), {B, 1), {A, 2), {B, 2), (1, 1). 

First, we consider about eq. ( ^.281) . For I = A,B, it is trivial from the definition of Ai 
and V 



M- 
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6 

Jj, 


6 

V 

6J 


6 


Jj 


6 

V 

SJ 


6 


Jj 



B 



f 5^ 



-de 



B 



For / = 1, using 



(^w)i 

.5J 5J L 



'10c 



5 



+ 



5 



+ 



5 



.SJAiC) 5JBicj)J 5Ji(C,fT) 



5JAiO '5Ji(C,a) 5JBia) "5Ji(C,a) 



we have 



, SJJ \ 6JJJ, 

= -2((mI] dJul] +(mI] oJmI' 



-2V10c 



6 



6 



-dr 



6 



6 



6 



-dr 



6 



+ 



6 



.SJa{C) 5Ji(C,a) 6JB{a) 5Ji(C,a) 



On the other hand, 



\\6J 6Jj^ 



6j)b. 



— V— ^ 



(E.l) 
(E.2) 



(E.3) 
(E.4) 



(E.5) 



which is nothing but the r.h.s. of ( Pi).5|) . Similarly, we can show the validity of the formula 
for / = 2 by noticing the following identities 

2 ^ (5 

^ r 7 /A x^o^c(Ci, C2)-Y— — - 

j=idJA{Cj) dJi{C,ai) 



6 



6Ja{0 ^5Ji(C,ai) 



SJAiCr 



5Ji(C,ai)' 



^i5Ja{Cj) ^Ja(C) V ^ 5Ja{C), 



. (E.6) 
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Next, we consider eq. (|5.29| ). Note that takes the upper-triangular form: A4ij = 
for I <J. For (/, J) = (A, A), 



J2MakMal(a^] (C;C'), (e.7) 

K,L \ dJ J K,L 



because J^al = ^a,l- For (/, J) = (A, 1), 



5J , 



^9c.%z^.(Ci,C')-/— - d^,d^,D,{C,X')-Y^—. (E.8) 



On the other hand, 



K,L \ (>J / K,L 



= y^cUX] iC-Xi)+U^) (C';Ci,^i), (e.9) 

which is nothing but the r.h.s. of ( |E.8| ). 

Similarly, for (/, J) = (^4,2), (1, 1), checking the formula is straightforward by utilizing 
the identities such as 



SJa{w) 

= %D,(C', Ci)DUz, C2)-J— + d^^D^iC, C2)DUz, Ci)-j^, (E.IO) 

dJA{w) oJa[w) 

D,{Ci,C[)DUCiX[)Ddz,w)-^ = d^A'D.iCiXi)-^,- (E.ll) 

6Ja{C) oJa{z) 

The validity for (/, J) = (5, 1), (5, 2) is obvious from the symmetry of M. with respect to 
A^ B. 
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